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INTRODUCTION. 



On peut définir la fonction gamma, soit, d'après les procédés de l'an- 
cienne Analyse, au moyen d'une expression déterminée, soit, confor- 
mément aux idées modernes sur la théorie des fonctions, en partant de 
certaines équations fonctionnelles. Si l'on fait abstraction de cette der- 
nière méthode qui n'a donné naissance qu'à de rares travaux ( x ), très 
importants d'ailleurs, on se trouve en présence de deux définitions, dues 
l'une et l'autre à Euler. La première, fondée sur la considération de la 
limite d'un produit, a été préconisée par Gauss( 2 ) et Liouville ( 3 ). La 
seconde, où T(x) est l'expression d'une intégrale définie, a été adoptée 
successivement par Euler, Legendre et presque tous les analystes. On 
doit chercher, sans doute, la raison de cette préférence exclusive dans 
les nombreux rapports qui relient l'étude de T(x) à celle des intégrales 
définies. Cependant la définition choisie par Gauss, non seulement pos- 
sède l'avantage d'une plus grande généralité, puisque la variable n'y 
est astreinte qu'à la seule condition restrictive de ne pas être égale à un 
entier négatif, mais encore elle révèle immédiatement la nature même 
de cette transcendante et permet d'établir toutes ses propriétés d'une 
manière plus concise, plus rigoureuse et aussi plus naturelle ; au lieu de 
reposer sur une suite d'artifices, parfois compliqués, les démonstrations 
se développent avec une remarquable uniformité. 

Les considérations précédentes nous sont inspirées par un Mémoire 
critique d'un mathématicien allemand de grand mérite, le professeur 
Pringsheim, de Munich ( 4 ). Après avoir comparé les avantages et les 



( ) Les principaux sont ceux de Weierstrass, Hankel, Prym, Mellin. 

( 2 ) Briefwechsel zwischen Gauss und Bessel, p. 152. 

( 3 ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. XXXV, 1852, 
p. 320. 

( 4 ) Zur Théorie der Gamma-Functionen (Mathematische Annalen, t. XXXI, 1888, p. 455- 
481). 



VI 



inconvénients des deux définitions dans le sens que nous venons d'indi- 
quer, ce savant auteur arrive à la conclusion, qu'au point de vue didac- 
tique, il est préférable, à l'exemple de la plupart des analystes, de définir 
T(x) au moyen d'une intégrale, puis d'en déduire que cette fonction est 
la limite d'un produit, et, ce résultat une fois obtenu, de l'utiliser pour 
l'exposé de la théorie. Cette manière de voir s'appuie assurément sur 
des arguments sérieux, et pourtant, malgré leur valeur incontestable, 
nous croyons devoir formuler un avis absolument opposé. Nous réfé- 
rant à l'opinion de Gauss et de Liouville, nous estimons, au contraire, 
que l'expression de T(x) sous forme de limite d'un produit s'impose dès 
le début et qu'il vaut mieux ne faire intervenir les intégrales définies 
qu'en dernier lieu, lorsque leur introduction devient indispensable. Il 
nous paraît, en effet, peu rationnel de recourir tout d'abord à un ordre 
d'idées que l'on abandonnera presque aussitôt pour le reprendre ensuite. 
D'ailleurs, n'est-il pas plus logique de choisir celle des deux définitions 
qui nécessite les notions les moins élevées et permet, par là même, de 
réaliser ce maximum de simplicité qui doit être le but primordial de 
toute investigation mathématique ? Dans ce travail essentiellement syn- 
thétique, nous nous sommes proposés de justifier cette manière de voir, 
en considérant spécialement le cas où la variable est réelle, sans ce- 
pendant négliger les idées récentes qui ont si profondément modifié la 
théorie de la fonction gamma. Nous n'avons pas cherché à découvrir des 
résultats nouveaux, mais nous avons apporté un soin scrupuleux dans le 
choix des méthodes ; en outre, la forme sous laquelle les raisonnements 
sont présentés nous est toujours personnelle ; souvent même, nous avons 
réussi à introduire des améliorations notables et nous le devons surtout 
à l'emploi presque constant des séries dont l'usage est toujours pra- 
tique et sûr. Enfin, quelque peu enclin à l'érudition mathématique, goût 
jadis si peu cultivé et qui maintenant se répand de jour en jour davan- 
tage, nous avons tenu, sans tomber dans une documentation excessive 
et fastidieuse, à donner incidemment de nombreux renseignements his- 
toriques et bibliographiques. Quelques-uns serviront à préciser certains 
faits ; tous, sans doute, contribueront à augmenter l'intérêt de la théorie. 



Maurice Godefroy. 



LA FONCTION GAMMA; 

THÉORIE, HISTOIRE, BIBLIOGRAPHIE. 



I. 

HISTORIQUE. 



L'introduction de la fonction gamma dans l'Analyse est due à Euler. 
Cependant, à propos de l'origine de cette transcendante, il serait injuste 
de ne point rappeler les résultats obtenus par Wallis et Stirling. Le pre- 
mier, par la découverte de sa fameuse formule ( x ), montra quel parti on 
pouvait tirer des quadratures pour l'évaluation de certains nombres ; le 
second, en cherchant à déterminer le logarithme du produit des n pre- 
miers entiers, parvint à une série des plus singulières dont l'étude se 
rattache étroitement à celle de la fonction gamma ( 2 ). 

Les premiers travaux d'Euler sur cette matière remontent à l'an- 
née 1729. Goldbach et Daniel Bernoulli s'occupaient alors de l'inter- 
polation de la série 

l + l-2 + l-2-3 + --- + l-2---n + --- 

considérée auparavant par Wallis. Euler, mis au courant de ce problème 
par Bernoulli, échangea avec Goldbach une correspondance ( 3 ) qui devait 



( ) Arithmetica infinitorum (Johannis Wallis. . . opéra mathematica, t. 1, p. 469). 

( )Methodus differentialis : sive Tractatus de summatione et interpolatione serierum infi- 
nitarum. Londini, 1730, in-4°. «Cet ouvrage, d'un génie analytique supérieur, offrait 
la solution de problèmes qui semblaient dépasser les méthodes de la science à cette 
époque». (Binet). 

( 3 ) Voir P. -H. Fuss, Correspondance mathématique et physique de quelques célèbres géomètres 
du xvin eme siècle. . . , 1. 1. 



se prolonger jusqu'à la mort de ce dernier survenue en 1764. Dans une 
première lettre du 13 octobre 1729, il observe que le terme général de la 
série de Wallis est la limite de l'expression 

J ' 2 ' 3 '" n -(n + l) m 



(1 + m){2 + m) • " (n + m) 



lorsque n croît indéfiniment ( 1 ). Cette question fait également l'objet 
d'une lettre du 8 janvier 1780 ( 2 ) ; il y annonce que le terme général de la 
même série peut être représenté par l'intégrale 



I 



dx(-lx) n . 



Ainsi, dès le début de ses recherches sur les fonctions qu'il qualifia plus 
tard d'inexplicables ( 3 ), Euler avait considéré le produit retrouvé plus tard 
par Gauss, mais il ne reconnut point quels avantages il offrait, et dé- 
sormais, dans les nombreux travaux qu'il publia sur ce sujet, il s'en tint 
uniquement aux intégrales définies comme mode de représentation. 

Le Mémoire de Gauss sur la série hypergéométrique ( 4 ) marque un 
nouveau progrès dans l'histoire de la fonction gamma. C'est un exposé 



( ) «Hujus seriei 1, 2, 6, 24, 120, etc., quam a Te multum tractatam esse vidi, hune 
inveni terminum generalem 



i ™ jl—mnm nl—mAm Al—mzm 

etc. 



l + m 2 + m 3 + m 4 + m 

ex infinito factorum numéro constantem, qui terminum ordine m mam exprimit. Is qui- 
dem in nullo casu abrumpitur, et aeque si m est numerus integer, tantum ad verum 
magis magisque accedit, ac si m fuerit fractus. Sed tamen per eum admodum prope 
quemque terminum invenire licet, idque eo facilius, quo minus assuma tur m. Si au- 
tem aliquot solum, uti visum sit, factoribus, termine» generali commodior induci potest 

1 • 2 

forma : ut si duobus prioribus factoribus contenti esse velimus, habebitur -3 m 

r (1 + m)(2 + m) 

pro termino ordine m. Sin autem generaliter n factures capiantur, sequentibus reliquis 

neglectis, erit terminus generalis — r — An + l) m , qui, quo major accipitur 

numerus n, eo propius ad verum accedet.» P.-H. Fuss, Correspondance mathématique et 
physique. . . , 1. 1, p. 3-4. 

( 2 ) P.-H. Fuss, Correspondance mathématique et physique. . . , 1. 1, p. 11-18. 

( 3 ) Le Chapitre XVI de la seconde partie des Institutiones Calculi differentialis, p. 769- 
807, est tout entier consacré aux fonctions inexplicables. «Il n'est pas facile d'expliquer 
comment Euler a pu appliquer une telle épithète (inexplicabïlis) à une grandeur dont 
il connaissait, depuis vingt-cinq ans, l'expression analytique sous forme continue». 
(Binet). 

( 4 ) Disquisitiones générales circa seriem infinitam 1 + — -x + ■ ■ ■ (Cari Friedrich Gauss 

Werke, t. III, p. 123-162). Gauss présenta son Mémoire à la Société royale des Sciences 
de Gôttingue le 30 janvier 1812. 



très sobre des propriétés les plus essentielles de cette transcendante re- 
gardée comme la limite, pour k - oo, du produit 

1 • 2 • 3 • • • k 

"^ = (2 + l)(2 + 2)(2 + 3)...(z + k)* Z(1) - 

Gauss, dans une lettre adressée à Bessel en date du 21 novembre 1811 ( 2 ), 
a expliqué les motifs qui l'ont conduit à choisir cette définition ( 3 ). Ce- 
pendant, et cela nous paraît assez étrange, il abandonna plus tard sa 
manière de voir primitive, car l'exemplaire imprimé de son Mémoire, 
dont il avait l'habitude de se servir, porte une annotation manuscrite où 
une opinion toute différente est exprimée ( 4 ). 

C'est à Legendre que l'on doit la première étude complète et appro- 
fondie des intégrales eulériennes. Dans ses Exercices de Calcul intégral ( 5 ), 
puis dans son Traité des fonctions elliptiques ( 6 ), se trouvent rassemblés 
tous les résultats acquis à cette époque dans cette partie de l'Analyse. 
Les démonstrations qui en sont données ont presque toujours, sinon 
le mérite de la rigueur, du moins celui de l'originalité, et souvent des 
développements nouveaux les accompagnent. 

Cependant les travaux consciencieux de Legendre étaient loin d'être 
définitifs; depuis lors l'attention des analystes s'est maintes fois fixée 
sur cette région restreinte, mais si attrayante, de l'immense domaine 
de la Science mathématique. Aussi, la simple énumération de tous les 
Mémoires écrits sur ce sujet depuis le commencement du siècle dernier 
suffirait à étendre démesurément cette courte notice. La plupart ont pour 
auteurs des savants de premier ordre. Nous rappellerons, entre beaucoup 
d'autres, les noms de Binet, Lejeune Dirichlet, Liouville, Cauchy Raabe, 



( ) Cari Friedrich Gauss Werke, t. III, p. 144-146. 

( ) Briefwechsel zwischen Gauss und Bessel, p. 152. 

( ) Voir aussi l'analyse que Gauss a faite de son Mémoire dans les Gôttingische gelehrte 
Anzeigen (Cari Friedrich Gauss Werke, t. III, p. 200-201). 

( ) Cette annotation, rapporte Schering, est relative au paragraphe 28 ; elle est ainsi 
conçue : «Die beste Définition von II ist dass 

Tlm = e (<* +1 V" dx». 



-f 



Voir : Cari Friedrich Gauss Werke, t. III, p. 230. 

( ) Exercices de Calcul intégral sur divers ordres de transcendantes et sur les quadratures, 
Paris, V e Courcier, 1811-1819, 3 vol. in-4°. Cet ouvrage est devenu rare. 

( 6 ) Traité des fonctions elliptiques et des intégrales eulériennes avec des Tables pour en faciliter 
le calcul numérique, Paris, Huzard-Courcier, 1825-1828, 3 vol. in-4°. Le tome II contient 
un Traité des intégrales eulériennes qui est la reproduction améliorée des parties II et IV 
des Exercices. 



Gudermann, Kummer, Weierstrass, Hermite, Prym, Appell, Lerch, Mel- 
lin. A cet égard, on trouvera des renseignements bibliographiques très 
détaillés dans l'excellente monographie de Brunel( x ). 

Si la fonction gamma a été l'objet de multiples travaux, il faut en at- 
tribuer surtout la raison à ce que l'on peut suivre graduellement et sans 
effort, dans l'exposé de ses propriétés, tous les progrès de l'Analyse mo- 
derne. En effet, type par excellence de la fonction uniforme, elle se prête, 
d'une manière particulièrement heureuse, à l'application des théorèmes 
fondamentaux de la Théorie générale des fonctions dont elle constitue 
peut-être l'illustration la plus simple et la plus saisissante. Nous avons vu 
que Gauss, se plaçant à un tout autre point de vue qu'Euler et Legendre, 
la définit comme la limite d'un produit convergent. La voie était ainsi 
ouverte à Weierstrass ; ce dernier fut effectivement amené à la décou- 
verte de son théorème sur les facteurs primaires par la considération de 
la fonction Yl(z) de Gauss, comme il le dit lui-même ( 2 ), dans son célèbre 
Mémoire Zur Théorie der eindeutigen analytischen Functionen ( 3 ) qui parut 
en 1876. D'ailleurs dans un travail bien antérieur ( 4 ), où il s'était proposé 
de traiter avec clarté et rigueur la question si complexe des facultés analy- 
tiques, Weierstrass avait reconnu le caractère propre de l'inverse de T(x) 
qu'il appela factorielle de x. 

Une autre preuve de l'influence constante que les perfectionnements 



(^Monographie de la fonction gamma (Mémoires de la Société des Sciences physiques et 
naturelles de Bordeaux, 3 e série, t. III, 1886, p. 1-184). Le travail de Brunel se recommande 
surtout par l'abondance des données bibliographiques, mais au point de vue didac- 
tique il est un peu succinct. Depuis sa publication, il a paru plusieurs monographies 
consacrées à l'étude de la fonction gamma; nous nous plaisons à les citer, car leur 
lecture nous a été grandement profitable. Ce sont, d'après leur ordre d'apparition : la 
dissertation inaugurale d'Arvid Lindhagen, Studier ôfver gamma-funktionen och nâgra 
beslagtade transcender, Upsala, 1887, in-4° ; les divers articles de J. L. W. V. Jensen insérés 
dans le Nyt Tidsskrift for Mathematïk, t. II B, 1891, p. 33-56, 57-72, 83-85, sous le titre 
de Gammafunktionens Theori i elementxr Fremstilling ; et enfin la Théorie funkce gamma 
de Mathias Lerch écrite en langue tchèque (Vëstnîk ceské Akademie cisafe Frantiska Jo- 
sefa pro vëdy, slovesnost a umënî v Praze, t. II, 1893, et t. VIII, 1900). Il faut mentionner 
aussi un Mémoire récent de Barnes, The theory ofthe gamma function (The Messenger of 
Mathematics, vol. XXIX, 1899, p. 64-128) qui est conçu d'après un plan assez spécial. 

( ) «Aber eben dièse Function weist auf den Weg him, der zum Ziele fùhrt» (Mathe- 
matische Werke, t. II, p. 91). 

( ) Mathematische Werke, t. II, p. 77-124. Une traduction de ce Mémoire par Picard 
a été insérée dans les Annales scientifiques de l'École normale supérieure, 2 e série, t. VIII, 
1879, p. 111-150. 

( ) Ùber die Théorie der analytischen Facultdten (Mathematische Werke, t. I, p. 153-221). 
Ce travail fut publié en 1856 dans le tome 51 du Journal fur die reine und angewandte 
Mathematïk (p. 1-60). La théorie des facultés analytiques avait été auparavant étudiée 
par Vandermonde, Kramp, Bessel, Crelle, Eytelwein, Oettinger, Ohm, Raabe, mais elle 
renfermait encore de graves lacunes. 



de l'Analyse ont exercé sur le développement de la théorie de la fonction 
gamma nous est encore fournie par les fonctions P(x) et Q(x). Ces fonc- 
tions introduites par de Gasparis ( x ), mais dont la signification véritable 
a été mise en relief par Prym ( 2 ), se présentent comme une conséquence 
fort remarquable du théorème de Mittag-Leffler. Un champ, jusqu'alors 
inexploré, était par là même livré à l'activité des chercheurs, et, parmi les 
plus beaux résultats obtenus, il importe de mentionner ceux qui ont été 
trouvés par Hermite ( 3 ) en appliquant le théorème de Mittag-Leffler 
à certaines combinaisons de fonctions gamma. Aussi bien, Hermite a 
publié de nombreux Mémoires sur les fonctions eulériennes, pour les- 
quelles il semble avoir éprouvé une prédilection marquée ; une part pré- 
pondérante lui est ainsi assurée dans l'avancement de cette théorie qui 
conservera toujours l'empreinte ineffaçable de son puissant génie. 

Tel est, en quelques mots, l'historique de la fonction gamma ; bien que 
brièvement esquissé, il suffit, sans qu'il soit besoin d'insister davantage, 
à montrer combien privilégiée est la place que cette transcendante tient 
en Analyse et quelle importance on est en droit d'attacher à son étude. 



1 î 

( ) Sul calcolo del valore délia funzione £ — — (Società reale di Napoli. Rendiconto dell'Acca- 

demia délie Scienze fisiche e matematiche , t. VI, 1867, p. 276-388. — Giornale di Matematiche, 
t. VI, 1868, p. 16-23). Si, comme le remarque Hermite, l'intérêt de la décomposition de 
T(x) en fonctions P(x) et Q(x) demeura caché à de Gasparis, c'est que les principes de la 
théorie des fonctions étaient alors peu connus. 

( ) Zur Théorie der Gammafunction (Journal fiir die reine und angewandte Mathematik, 
t. LXXXII, 1877, p. 165-172). 

( ) Sur une application du théorème de M. Mittag-Leffler, dans la théorie des fonctions 
(Journal fur die reine und angewandte Mathematik, t. XCII, 1882, p. 145-151). Voir aussi : 
Faculté des Sciences de Paris. — Cours de M. Hermite rédigé en 1882 par M. Andoyer, 3 e éd., 
p. 134-138. Le Mémoire de Mellin, Zur Théorie der Gammafunction (Acta Mathematica, 
t. VIII, 1886, p. 37-80), mérite d'être placé à côté des travaux de Hermite. 



IL 



ETUDE GENERALE DE LA FONCTION GAMMA. 



Soit 



n(x) = ri 



1-2-3 



x(x + 1) • •• (x + ri)' 

il est facile d'établir que, pour toute valeur de x non égale à zéro ou à un 
entier négatif, cette expression tend vers une limite quand n augmente 
indéfiniment. En effet, de l'identité 



2-3- 



n - 



1-2- ■■{n-1)' 



on tire immédiatement 



-"-"♦ïiKr 



i+ 



n - 1, 
de sorte que le produit n(x) peut être mis sous la forme 

1 \X 



un * K)vr 

n(x) = ; — 



1 + 



n-1 



1 + ïA 1 + 2 



1 + 



si l'on prend alors les logarithmes, en adoptant, pour éviter les valeurs 
multiples, la détermination de la fonction logarithmique qui s'annule 
lorsque la variable est égale à l'unité, on trouve que le logarithme du 
produit xll(x) est égal à 



p=n r -, p=n r -, 



expression dans laquelle la seconde somme se déduit de la première en 
faisant dans celle-ci x = 1. Mais, dès que n dépasse le plus grand entier 
contenu dans le module de x, on a 



£-jog( 

n \ 



! + -)=-(" 

ni 2 \n 



1 (x 

4 In 



par suite, si r est un entier positif arbitrairement grand, lorsqu'on rem- 
place dans ce développement x par r, et tous les coefficients par l'unité, 
on obtient l'inégalité 

v 2 



X ! (t X 

--log 1 + - 
n \ n 



< 



n(n - r)' 
qui est vérifiée à partir de n - r + 1 si le module de x reste inférieur à r. 



Ainsi la série 



i-logl + i 



+ 



x . /, , x 



+ 



est absolument et uniformément convergente pour toute valeur de x inté- 
rieure au cercle de rayon r et non égale à un entier négatif. En particulier, 
pour x - 1, la série positive 



i-^KH^K) 



+ 



est convergente ; sa somme, que nous représenterons par p, est connue 
sous le nom de constante d'Euler{ 1 ). 

Il résulte de ce qui précède que le produit n(x), pour toute valeur de x 
autre que zéro ou qu'un entier négatif, a une limite lorsque n devient 
infini. Cette limite est une fonction transcendante de x à laquelle on 
donne le nom de fonction gamma ( 2 ) ; on la désigne, en effet, par le symbole 
T(x) employé pour la première fois par Legendre ( 3 ), et unanimement 



( 1 ) Commentant Academix Scientiarum imperialis petropolitanse, t. VII, 1734-1735, p. 156- 
157. 

( ) Il est facile d'imaginer des fonctions plus générales que T(x) définies au moyen 
d'une expression analogue à n(x). Voir à ce sujet une Note de Cesàro dans les Nouvelles 
Annales de Mathématiques, 3 e série, t. IX, 1890, p. 356-359. On peut également généraliser 
la fonction log T{x) ; cette question a été traitée par Berger, Ôfversigt afKongl. Vetenskaps- 
Akademiens Fôrhandlingar, t. XXXVIII, 1881, n° 9, p. 13-30. 

( ) Mémoires de la classe des sciences mathématiques et physiques de l'Institut de France, 
année 1809, p. 477. Gauss, comme nous l'avons déjà dit (p. 3), s'est servi des notations 
ll(/c, z) et ll(z) pour représenter le produit 

l-2-3---k 

(z + l)(z + 2)(z + 3) ■ ■ ■ (z + k) 

et sa limite lorsque k devient infini. Ainsi la fonction Il(z) de Gauss ne diffère pas de la 
fonction T(z + 1) de Legendre. 



conservé depuis. 

La constante d'Euler qui, comme on vient de le voir, est la limite, pour 
n - oo, de l'expression 

p n = s n -logn, 

où s n représente la somme des n premiers termes de la série harmonique, 
joue, dans la théorie de la fonction gamma, un rôle analogue à celui du 
nombre n dans la théorie des fonctions circulaires^); sa valeur a été 
calculée avec 59 décimales exactes par Shanks ( 2 ) ; en se bornant à vingt 
décimales, on a 

p = 0.57721566490153286060 .... 

Weierstrass ( 3 ), dans son Mémoire sur les facultés analytiques, a donné 
le nom de factorielle de x à l'inverse de T(x), et il a mis cette fonction, qu'il 
représente par Fc(x), sous une forme que nous allons faire connaître, car 
nous aurons l'occasion de l'utiliser. 

Si l'on considère l'égalité 

1 1 x(x + !)••• (x + n) 





U(x) n x 1 • 2 • • • n 


/ 


on en déduit 




U(x) \ n) 


[¥©1 


~x + 2 (2\ x ~ 
2 \3/ 




x + n - 1 (n - l\ x 
n -1 \ n I 



et, pour n - oo, on obtient 

M=oo 

Fcw-xn(i+; M „ +1 

M = l 

telle est l'expression d'où Weierstrass a tiré toutes les propriétés de Fc(x). 



Formule de Weierstrass. — On a 



U(x) 



1 x(x + l)(x + 2) • • • (x + n) 
n x 1 • 2 • • • n 



( ) On peut consulter, pour l'histoire de ce nombre célèbre, la monographie publiée 
par Glaisher dans The Messenger of Mathematics, vol. I, 1872, p. 25-30, et vol. II, 1873, 
p. 64. Nous avons donné, dans notre Théorie élémentaire des séries, une étude sommaire 
de la constante d'Euler. 

( ) Proceedings ofthe Royal Society ofLondon, 1869-1870, vol. XVIII, p. 49. 

( ) Mathematische Werke, 1. 1, p. 161. 



ou bien, s n désignant la somme des n premiers termes de la série harmo- 
nique, 



Yl(x) 



1 e n 



i + fW-f 



i + ^-ï 



1 x \ -- 

1 + -)e n 

n 



Soit 



gSft 

A n - — ; 
n 



en prenant les logarithmes, on trouve 

logÀ„ = s„-logn; 

quand n croît indéfiniment la limite du second membre de cette égalité 
est la constante d'Euler p, par suite 



T(x) 



xe 



px 



i + f u-f 



l^fi 



1 + - ]e n 

n 



La relation précédente, qui peut servir à définir la fonction gamma, 
donne la décomposition de son inverse en facteurs primaires. Elle a été 
signalée par Weierstrass ( x ) dans son Mémoire sur la théorie des fonctions 
uniformes. 

Si l'on se reporte à la formule 



sm7ix 



TCX 



n=oo n=oo 

1 + - )e n 

ni 
n=\ n=l 



n 



n 



X\ * 
-)en, 

n) 



on constate que la fonction T(x) est formée avec la moitié des facteurs qui 
constituent la fonction simplement périodique sin nx. Ce fait a conduit 
à imaginer des fonctions formées avec la moitié, ou même le quart, 
des facteurs qui constituent certaines fonctions doublement périodiques. 
Telles sont les fonctions de Heine ( 2 ) et la fonction gamma double de 
Barnes ( 3 ). A ces fonctions se rattachent aussi les fonctions introduites 



( ) Mathematische Werke, t. II, p. 91. Cependant cette formule avait été déjà rencontrée 
par Newman (The Cambridge and Dublin mathematical Journal, vol. III, 1848, p. 59). 

( 2 ) Journal fur die reine und angewandte Mathematik, t. XXXIV, 1847, p. 309-315. — 
Handbuch der Kugelfunctionen, 2 e éd., 1. 1, p. 109-115. 

( 3 ) Proceedings ofthe London mathematical Society, vol. XXXI, 1899, p. 358-381. — The 
quarterly Journal ofpure and applied Mathematics, vol. XXXI, 1900, p. 264-314. 
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par Appell ( x ), par Pincherle ( 2 ) et par Alexejewsky ( 3 ). 

La formule de Weierstrass montre que la fonction T(x) est méro- 
morphe pour toute valeur finie de x. Elle admet, en effet, comme seuls 
points singuliers, des pôles simples correspondant aux valeurs suivantes 
de la variable : 

x - 0, x = -1, x - -2, . . . , x - -n, ... r 

le point x - oo étant un point singulier essentiel. 

Théorème. — L'inverse de la fonction gamma est développable en une série 
entière de rayon de convergence infini. 

Il résulte immédiatement de la décomposition de l'inverse de la fonc- 
tion gamma en facteurs primaires que cette fonction est développable 
en une série entière convergente dans tout le plan. Cette proposition, 
énoncée pour la première fois par Weierstrass ( 4 ), peut d'ailleurs être 
démontrée directement de la manière suivante : 

Si l'on pose 

(x \ _ x 
1 + -je n, 

on en déduit 

l_x*_ 2 x 3 i\p _Pzl %P 

Un ~ l-2 n 2 + l-2-3n3 '" ( } l-2---pnP '"' 

Soit alors r le module de x ; si l'on considère la série positive 
\ r 2 2 r 3 p-1 rP 

+ - — ^~z -=■ + •••+ — — + • • • , 



" 1 • 2 n 2 1-2-3 n 3 1 ■ 2 • • • p nP 

ses coefficients étant moindres que l'unité, dès que n atteint l'entier m 
immédiatement supérieur à r, on a 



a n < 



r 2 



n(n - r)' 
d'où il suit que la série positive de terme général a n est convergente. 



(*) Mathematische Annalen, t. XIX, 1882, p. 84-102. 

( ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. CVI, 1888, 
p. 265-268. 

( ) CooSmeHia h npoTOKOJiJibi XapLKOBCKaro MaTeMaTH^jecKaro OGmecTBa, 2 e sé- 
rie, t. I, 1890, p. 169-238. — Berichte tiber die Verhandlungen der kôniglich sachsischen 
Gesellschaft der Wissenschaften zu Leipzig. Mathematisch-physische Classe, t. XLVI, 1894, 
p. 268-275. 

( ) Ueber die Théorie der analytischen Facultàten (Mathematische Werke, 1. 1, p. 161). 
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Le produit infini de ternie général l+u n , d'après un théorème connu if), 
est donc développable en une série entière convergente pour toute va- 
leur de la variable; par conséquent, l'inverse de T(x) est une fonction 

transcendante entière. 

I 

La fonction — -. pour x - 0, se réduit à l'unité; autrement dit, le 

r(i + x) 1 

premier terme de son développement n'est pas nul ; on en conclut que 
la fonction T(l + x) est développable suivant une série entière de rayon 
de convergence non nul ( 2 ) ; la valeur de ce rayon de convergence est 
d'ailleurs égale à l'unité, puisque le point singulier le plus rapproché de 
l'origine est le point x - -1. 

Les coefficients des développements de la fonction T(l + x) et de son 
inverse seront calculés ultérieurement. 

Relation fonctionnelle. — On a 

, __ x+ \ l-2-3---n _ nx 

( ) = n (x + l)(x + 2)---(x + n + l) = n + x + 1 W ' 

d'où, pour n - oo, 

r(x + l)=xr(x); 

c'est cette relation que l'on nomme relation fonctionnelle. Elle ne caractérise 
pas d'ailleurs la fonction gamma, car elle est vérifiée par toute fonction 
de la forme 

y - r(x)0(x), 

la fonction 0(x) satisfaisant elle-même à l'égalité 

0(x + 1) - 0(x), 

c'est-à-dire admettant une période égale à l'unité, ce qui est le cas de 
toutes les fonctions développables en séries trigonométriques conver- 
gentes telles que 

0(x) = «g + a i cos 27zx + «2 cos 47TX h + b\ sin 2nx + 1)2 sin Aux H . 

Si dans la relation fonctionnelle 

T(x + 1) = xT(x), 

on remplace successivement x par x + 1, x + 2, . . . , x + m - 1, on trouve 

T(x + 2) - (x + 1) T(x + 1), 



T(x + m) = (x + m - 1) T(x + m - 1), 



() Jules Tannery, Introduction à la théorie des fonctions d'une variable, p. 206-210. 
( ) Maurice Godefroy, Théorie élémentaire des séries, p. 96. 
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d'où, m étant un entier positif, 

T(x + m) = x(x + 1) • • • (x + m - 1) T(x). 

Soit x une valeur réelle de la variable extérieure à l'intervalle (0, 1) ; en 
désignant par m l'entier positif immédiatement inférieur à x ou supérieur 
à -x suivant que x est positif ou négatif, la formule précédente permet 
de ramener le calcul de T(x) à la détermination de cette fonction pour 
une valeur de x intérieure à l'intervalle (0, 1). En effet, si x est positif, on 
a 

T(x) - (x - m)(x - m + 1) • • • (x - 1) T(x - m), 

et si x est négatif, on peut poser 



T(x) = 



x(x + 1) • • • (x + m - 1) 



T(x + m). 



Lorsque la variable x est imaginaire, ces deux relations donnent le 
moyen de réduire le calcul de T(x), pour une valeur quelconque de x, 
à l'évaluation de T(x) pour une valeur de x dont la partie réelle soit 
intérieure à l'intervalle (0, 1). 

Si la variable est un entier m + 1 on obtient simplement 

r(w + l) = 1-2- ■■m. 



Applications. — I. Module de T(a + §ï). — On vient de voir que si la 
variable est imaginaire, on peut toujours ramener sa partie réelle à être 
comprise entre et 1. Soit donc 

x - a + f>i, 
le nombre a satisfaisant aux conditions 

S a ^ 1; 
on a 

w ^— t = Xû^Tf f(l + a -) 2 + tfe-i, 
\T(a + Bi)\ y H 11 l ni n 2 

n=l L J 



mais 



n 



1 a ? ^ 
1 + - +^ 

n) n 1 



1 

2 _a 
e n 



n=co 



n=\ 



n=co 



0KHD 



M=l L 



1 + 



^ 



(a + n) 2 



13 



d'ailleurs 



r(l + a) 



n=co 

= e a P n( i + ~^ - "' 

n=l 



n 



et, quant au produit 



n=oo 

n 



1 + 



F 



il est compris entre les deux produits 



{a + n) 2 



n 

M = l 



1 + 



F 



(n + l) 2 



n 






qui correspondent aux valeurs extrêmes de a. Ces deux produits ont 
respectivement pour limites 



1 e$> n - e'P 71 



1 + f lf>n ' 



on peut donc poser 



|r(a + j8i)| = A 



r(l + a) 



2JStz 



2jStt 



2 iefr-e-P 71 ' 



le coefficient À vérifiant la double inégalité 

1 < A < -^1+jS 2 . 

Cet élégant résultat est dû à Mathias Lerch ( x ). 

Lorsque a est nul, À se réduit à l'unité et l'on retrouve la formule 



TOI 

qui a été donnée par Stieltjes ( 2 ). 



2/37T 



fi y e pn_ e -pn 



IL Résidus de T(x). — Le calcul du résidu A n correspondant au pôle 
x - -n est bien simple. Dans le domaine de ce point, on peut poser 



T(x) = F(x + n) + 



x + n' 



(}) Théorie funkce gamma (Vëstntk ceske Akademie cisafe Frantiska Josefa pro vëdy, sloves- 
nost a umënî v Praze, t. II, 1893). 

( ) Journal de Mathématiques pures et appliquées, 4° série, t. V, 1889, p. 426. 
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en désignant par F(x + n) la somme d'une série procédant suivant les 
puissances entières successives du binôme x+n ; de la relation précédente 
on tire alors 

A n = (x + n) T(x) - {x + n) F(x + n). 

Le résidu A n est donc la limite, pour x - -n, du produit 

vw s T(x + n + 1) 

{X + n)T{x) = x(x + l)...(x + n-lY 

par suite 

A - ( - 1)M 



1-2---H 



1 T(x + n) 
Limite pour n - oo de l'expression — —^ — -. — On a 

r r n 1 1 ■ 2 ■ ■ ■ (n - 1) 

1 T(x + n) 1 x(x + 1) • • • (x + n - 1) 

■ 1 (x), 



n x i . 2 . . . ( n - 1) M * 1 • 2 • • • (n - 1) 

c'est-à-dire 

1 T(x + n) n T(x) 



n x 1 • 2 • • • (n - 1) ft + x n(x) ' 
la limite de l'expression 

1 T(x + n) 
n* l-2---(n-l) 

pour n - oo, est donc égale à l'unité. Ce résultat indiqué par Weier- 
strass ( 1 ), joint à la relation fonctionnelle, suffit évidemment à caractériser 
la fonction gamma ; car, le supposant vérifié, on en déduit sans peine, au 
moyen de la relation fonctionnelle, que T(x) est la limite du produit n(x). 
On a été ainsi amené, par analogie, à étudier des fonctions définies 
par la double condition de vérifier une équation semblable à la relation 
fonctionnelle de T(x), et de satisfaire à une autre relation quelconque. 
Nous nous bornerons à traiter le plus simple des problèmes de ce genre 
dont la solution est donnée par les fonctions de Prym. 

Applications. — La propriété de la fonction gamma qui vient d'être 
établie fournit parfois un procédé commode pour la discussion de cer- 
taines séries dont les termes sont formés avec des factorielles. Nous n'en 
citerons que deux exemples se rapportant à notre sujet : 



(!) Ùber die Théorie der analytischen Facultàten (Mathematische Werke, t. I, p. 166 et 
p. 192-194). 
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I. Étude d'une série de Stirling. — Soient 

a n = a(a + 1) • • • (a + n - 1), 
b n = b(b + l)---(b + n- 1), 

en convenant de regarder aq et b$ comme égaux à l'unité ; on a 



a n _ a n+l 
bn b n+ \ 



■■•> m a + n\ a n 

= ~ ! 1 " n~ = ( b ~ fl ) Z 

b + nl b n+1 



b n 



d'où 



a a 1 a 

b b 1 b 1 

b\ &2 b 2 



a n-l a n ,, v a n-\ 

h h~ = {b ~ } ~h~' 

v n _i n Vn 



par suite 



1 ( 

b -a \ 


a n \ _ a ai 
b n I b\ b 2 


+ 


a n-l 
b n 


Or 


a n 1 T(b) 


1 T(a + n) 




1 T(b + n) 


b n n b ~ a T(a) 


n a i-2---(n-l) 




n b 1 • 2 • • • (n - 1) 


si donc on pose 









a- a + fii, b = a' + fi'i, 



.b—a 



le module de n a étant égal à n a a , on voit que, n croissant indéfiniment, 
le rapport -^- augmente au delà de toute limite pour a' - a < et tend 

u n 

vers zéro pour a' - a > 0, tandis que pour a' - a = 0, il est indéterminé. 
On en conclut que le développement 



1 

- - + 



+ 



a(a + 1) 



b b(b + 1) b(b + ï)(b + 2) 



+ 



est valable en supposant positive la partie réelle de la différence b-a;de 
plus, cette condition étant remplie, il est absolument convergent, car 



n 



a'-a+l 



a n 



b n +i 



n T(b) 



n + b T(a) 



T(a + n) 



n a 1 • 2 • • • (n - 1) 



T(b + n) 



n b i.2--.(n-l) 
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égalité dont le second membre tend vers une limite pour n - oo. La série 
précédente est l'une de celles que Stirling( 1 ) a considérées. 

IL — Convergence de la série hypergéométrique. — Le terme gé- 
néral de la série hypergéométrique 

_ a(a + 1) • • • (a + n - 1) jS(jS + 1) • • • (j8 + n - 1) n 

1 • 2 • • • n y(y + 1) • • • (y + n - 1) 

est, comme on le constate aisément, le produit des deux expressions 



n a+p-y-l TQ 7 ) n 



et 



T(a + n) 1 r(j8 + n) 



n a l-2---(n-l) n P l-2---(n-l) 



1 r(y + n) 
n) 7 l-2---(n-l) 



cette dernière ayant pour limite l'unité lorsque n devient infini. D'ailleurs, 
si a, b, c sont les parties réelles respectives des paramètres a, f}, y, le mo- 
dule de l'exponentielle n a+ P~^ est n ; par suite, on peut poser 

\Un\ — "-n W |^| , 

en désignant par A n un coefficient qui tend vers une limite pour n = oo. 
Le rayon de convergence de la série est évidemment l'unité ; il reste à 
examiner la nature de la série pour \x\ = 1. Trois hypothèses sont alors à 
distinguer : 

1° a+b—c—1 > 0. — Les modules des termes augmentent indéfiniment. 

2° a + b — c — 1 - 0. — Les modules des termes tendent vers une limite 
non nulle. 

3°a + b-c-l<0. — Les modules des termes tendent vers zéro. 

La série des modules est convergente quand a, b, c vérifient l'inégalité 

a + b - c < 0, 
et seulement dans ce cas, car 

n \u n \ = A n (). 



(^Methodus differentialis : sive Tractatus de summatione et interpoïatione serierum infini- 
tarum, p. 12. 

( ) Ce résultat a été utilisé par M. Appell à propos de son théorème bien connu sur 
les séries entières {Comptes rendus hebdomadaires des séances de l'Académie des Sciences, 
t. LXXXVII, 1878, p. 692). 
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Fonctions de Prym. 

On peut se proposer de déterminer la forme générale des fonctions 
F(x) satisfaisant à la relation fonctionnelle 

F(x + 1) - xF(x) = A, 

où A désigne une constante, et telles que la limite, pour n - oo, de 
l'expression 

1 F(x + n) 

n* l-2---(n-l) 

soit égale à une autre constante B. 

La relation fonctionnelle, supposée satisfaite, donne successivement 

F(x) = - F(x + 1) -A -, 

F(x + 1) = F(x + 2) - A 1 



x + 1 x + 1' 

F(x + 2) = — *— F(x + 3) - A^— , 

x+2 x+2 



1 1 

F(x + n) = F(x + n + 1) - A- 



x + n x + n 

si l'on ajoute ces égalités multipliées respectivement par les rapports 



F(x + n + 1) A 
x(x + 1) • • • (x + n) 


1 1 1 

x x(x + 1) ' ' x(x + 1) • • • (x + n) 


F(x + n + 1) 


1 A + (x + n) F(x + n) 
— ri( y) 



1' x' x(x + l)' '"' x(x + 1)- • • (x + n + 1)' 
on trouve 

F(x) 

or 

x(x + 1) • •• (x + n) n n x l ■ 2- ■ ■ (n - 1) 

et, la limite de cette expression étant £>r(x), si l'on représente par S(x) la 
somme de la série convergente 

11 1 

x x(x + 1) x(x + l)(x + 2) 

on a 

F(x) - BT(x) - AS{x); 
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telle est la forme générale des fonctions F(x) ; la fonction gamma est la 
plus simple ; elle correspond au cas où A est nul et B égal à l'unité. 
Si maintenant, dans la relation fonctionnelle 

F(x + 1) = A + xF(x), 

on remplace x successivement par x - 1, x - 2, . . . , x - m, on en déduit 

F(x)=A + (x-l)F(x-l), 
F(x - 1) = A + (x - 2)F(x - 2), 
F(x - 2) = A + (x - 3)F(x - 3), 



F(x - m + 1) = A + (x - m) F(x - m); 

en multipliant ces égalités respectivement par les produits 

1, x-1, (x-l)(x-2), ..., (* - l)(;c - 2) • ■ • (* - m + 1), 

on obtient 

F(x) = A[l + (x-l) + (x- l)(x - 2) + • • • + (x - l)(x - 2) ■ ■ ■ (x - m + 1)] 
+ (x - l)(x - 2) • • • (x - m) F(x - m), 

d'où, pour x = m + 1, l'entier m étant positif, 

F(m + 1) = A[l ■ 2 ■ ■ ■ m + 2 ■ 3 • • • m + • • • 

+ (m- \)m + m + 1] + 1 • 2 • • • tn(B - Ae). 

La série S(x) est le produit de deux séries qu'il est facile de déterminer. 
Soit, en effet, 

1 Uq Cl-y a n 



x(x + l)(x + 2)- ■ -(x + n) x x + 1 x + n 

si l'on chasse les dénominateurs, puis que dans l'égalité obtenue on fasse 
x égal successivement à 0, -1, -2, . . . , -n, on trouve 



1 1 1 


1 1 


au ~l-2---n ai ~ ll-2-(n-l)' 


" z 1 • 2 1 • 2 • • • (n - 2) 


par suite 




1 1 1 


1 1 1 


x(x + 1) • • • (x + n) 1 • 2 • • • n x 


l-2---(n-l) 1 x + 1 


1 
+ 1 • 2 • • • (n - 


1 1 


■2) 1-2 x + 2 


+ (D" a 


1 



1 • 2- • -n x + n' 
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cette expression est précisément le terme général du produit des deux 

séries 

1 1 1 
e = 1 + 7 + 7-7; + ■■■+ 7-7, + ••• , 

11-2 1 -2---W 

et 

111 11 , _„ 1 1 
1 1- (-1)" 1 ■ 

x lx + 1 l-2x + 2 l-2---n x + n 

on désigne par P(x) la somme de cette dernière série, et par Q(x) la 
différence entre P(x) et T(x), de sorte que 

P(x) = -S(x), Q(x) = T(x)-P(x); 
e 

les fonctions P(x) et Q(x) sont les fonctions de Prym( 1 ). La fonction 
T(x) se présente alors comme la somme des fonctions P(x) et Q(x), et 
l'expression générale des fonctions F(x) devient 

F(x) = (B - Ae) P(x) + BQ(x). 

La fonction F(x) se réduit à l'une ou à l'autre des fonctions de Prym, 
suivant que l'on choisit pour les constantes arbitraires A et B les nombres 

A 1 = —, B 1 = 0, 
e 

ou 

A 2 = -, B 2 = 1; 

e 

il en résulte que les fonctions P(x) et Q(x) vérifient les relations fonction- 
nelles 

P(x + l) = xP(x)--, Q(x + l) = xQ(x) + -, 
e e 

et sont telles que les rapports 



P(x + n) 1 Q(x + n) 



n x l-2---(n-l)' n x l-2---(n-l)' 
tendent, pour n = 00, le premier vers zéro, le second vers l'unité. 



() Journal fur die reine und angewandte Mathematik, t. LXXXII, 1877, p. 165-172. L'ex- 
pression de Q(x) sous forme explicite est compliquée. Hermite {journal fur die reine und 
angewandte Mathematik, t. XC, 1881, p. 332-336) est parvenu à un développement où 
figurent plusieurs séries. Mellin (Acta mathematica, t. II, 1883, p. 231-232) a fait connaître 
un autre développement, un peu plus simple, où n'entre qu'une seule série. 
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Enfin on a 

eP(x) = -1 - (x - 1) - (x - l)(x - 2) (x - 1) • • • (x - m + 1) 

+ (x - l)(x - 2) • • • (x - m) e P(x - m), 

e Q(x) = 1 + (x - 1) + (x - l)(x - 2) + • • • + (x - 1) • • • (x - m + 1) 
+ (x - l)(x - 2) • • • (x - m) e Q(x - m), 

d'où, pour x = m + 1, l'entier m étant positif, 

e P(m + l) = -l-2---m-2-3---m m - 1 + 1 • 2 • • • me, 

e Q(m + 1) = l-2---m + 2-3---m + --- + m + l. 

La décomposition de la fonction gamma en fonctions de Prym ( a ), qui 
vient d'être obtenue par une voie très élémentaire, constitue aussi une 
application extrêmement simple du théorème de Mittag-Leffler. 

En effet, la fonction T(x) n'admettant comme points singuliers que les 
pôles simples 

JL Vj^ JL -L f JL ^f • • • f JL ïlf m m m f 

et le point singulier essentiel x = oo, si l'on considère la série 

11111 11 
P( X ) = _ _ _ 1 + _!^ - + (_!)" - — + ... 

x lx + 1 l-2x + 2 l-2---n x + n 

dont les termes sont les parties principales correspondant aux pôles, on 
voit qu'elle est uniformément convergente pour toute valeur de x autre 
que zéro ou qu'un entier négatif ; car, si r est un entier positif aussi grand 
que l'on veut, à partir de n - r + 1, l'inégalité 



(-lf 



1 • 2- • -n x + n 



< 



1 • 2- • -n n 



est vérifiée pour toute valeur de x située à l'intérieur du cercle de rayon r 
et non égale à zéro ou à un entier négatif. Ainsi la série précédente 
représente la partie méromorphe de T(x), tandis que la différence 

Q(x) = T(x)-P(x) 



( ) Cette décomposition, lorsque la fonction gamma est mise sous forme d'intégrale 
définie, revient au partage de l'intervalle d'intégration. Pour l'étude simultanée des 
fonctions T(x), P(x), Q(x), faite au seul point de vue des intégrales définies, on pourra 
se reporter à un Mémoire de Scheefer (Journal fur die reine und angewandte Mathematik, 
t. XCVII, 1884, p. 230-241). 
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est holomorphe dans tout le plan. On en conclut que P(l + x) est dévelop- 
pable suivant une série entière de rayon de convergence égal à l'unité, et 
que Q(x) est une fonction transcendante entière. 

On obtient également des résultats fort élégants en appliquant le 
théorème de Mittag-Leffler à des combinaisons de fonctions gamma. 
Hermite ( x ) a donné plusieurs exemples de ce genre. D'ailleurs, on peut 
aller beaucoup plus loin dans cet ordre d'idées, comme l'a prouvé Mel- 
lin ( 2 ), en étudiant très en détail la fonction 

F(x) = e mx ra (* + fl ) r ^ + k)r A (* + 



i a ' (x + a') r/ 3 ' (* + &') r A ' (* + /')' 

où m, a, b, . ..,/, a' ,b' ,...,V désignent des constantes quelconques, et 
a, j8, . . . , A, a' , §',..., A' des exposants entiers et positifs. Les propriétés 
des fonctions de Prym apparaissent alors comme des cas particuliers 
de celles de la fonction F(x), et l'emploi du théorème de Mittag-Leffler 
permet d'obtenir un grand nombre de transcendantes nouvelles qui ont 
des rapports étroits avec la fonction gamma. Enfin, en se plaçant tou- 
jours à ce même point de vue, on est naturellement amené à se poser 
la question suivante, qui nous a été indiquée par M. Appell, et dont le 
développement sortirait du cadre de ce travail : 

Étant donnée une fonction rationnelle r(x), on considère comme formant 
une même famille de fonctions toutes les fonctions définies par la relation 

F(x) = e mx — - K(tang nx), 

T a (x + a)TP(x + b)---T A (x + l) 

où R(tangnx) est une fonction rationnelle quelconque de tangnx et où les 
constantes m,a,b,...,l, a' , b' ,... , V , ainsi que les exposants a, fi,..., A, a' , 
§',..., A' sont déterminés de façon que 

F(x + 1) = r(x)F(x). 

Trouver, pour cette famille de fonctions, une formule de décomposition en élé- 
ments simples analogue aux formules connues de décomposition des fonctions 
doublement périodiques de deuxième et de troisième espèce qui admettent des 
multiplicateurs donnés. 



(^Journal fur die reine und angewandte Mathematik, t. XCII, 1882, p. 145-151 . Voir aussi : 
Faculté des Sciences de Paris. — Cours de M. Hermite rédigé en 1882 par M. Andoyer, 3 e éd., 
p. 134-138. 

( 2 ) Acta mathematica, t. VIII, 1886, p. 37-80. 
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Application. Fonction de Bourguet. — On vient de voir que, pour 
toute valeur de x non égale à zéro ou à un entier négatif, on a 

11 1 

eP(x) = - + — T7 + -, 777 ^ + '"' 

x x(x + 1) x(x + l)(x + 2) 



or, 



par suite 



T(x + n) - x(x + 1) • • • (x + n - 1) T(x), 



P(x) 1 1 1 

e f(x) ~ T(x + 1) + T(x + 2) + T(x + 3) + " ' ' 

Nous donnerons le nom de fonction de Bourguet à la somme de cette 
série, en mémoire du savant modeste et ingénieux, qui, le premier, l'a 
considérée ( x ). Si l'on représente par T(x) la fonction de Bourguet, en 
adoptant la notation de Weierstrass, on obtient 

T(x) = Fc(x + 1) + Fc(x + 2) + Fc(x + 3) + • • • , 

d'où 

T(x) - T(x + 1) - Fc(x + 1); 

d'autre part T(x + n) tend évidemment vers zéro pour n = oo ; ce résultat, 
joint à la relation fonctionnelle précédente, suffit à caractériser la fonction 
T(x), car des relations 

T(x) - T(x + 1) - Fc(x + 1), 
T(x + 1) - T(x + 2) - Fc(x + 2), 



T(x + n - 1) - T(x + n) = Fc(x + n), 
on déduit 

T(x) - T(x + n) = Fc(x + 1) + Fc(x + 2) + • • • + Fc(x + n), 

d'où, pour n = oo, 

n=oo 

T(x) - ^ Fc(x + n). 

M = l 

La fonction T(x) est holomorphe dans tout le plan, et par suite dé- 
veloppable suivant une série entière de rayon de convergence infini. 



( 1 )Acta mathematica, t. II, 1883, p. 280-287. Lindhagen a donné les tables de T(x) 
de -6 à 2, de centième en centième, et à trois décimales (Studier ôfver gamma- funktionen, 
p. 56-67). 
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Étude de l'équation P(x) - 0. — La formule 

e P(x) = - + — — + — — — + • • • 

x x(x + 1) x(x + l)(x + 2) 

montre que la fonction P(x) est positive pour toute valeur positive de x 
et décroît quand x augmente. Nous allons chercher maintenant le signe 
de P(x) dans un intervalle (-n, —n + 1) ; il suffit pour cela de déterminer 
le signe de eP(x - n), c'est-à-dire de la fraction 

x + (x - l)(x - 2) + • • • + (x - l)(x - 2) • • • (x - n + 1) + e P(x) 
(x - l)(x - 2) • • • (x - n) ' 

quand x varie de à 1. 

Si l'on considère d'abord la fraction 

„ N l + eP(x) 
eP(x-l)= x _l\ 

on constate immédiatement que la fonction P(x) est négative dans l'in- 
tervalle (-1,0). Au contraire, elle est positive dans l'intervalle (-2,-1), 

car 

x + eP(x) 

Si, dans le second membre de cette relation, on fait varier x successive- 
ment dans les intervalles (-1, 0), (-3, -2), (-5, -4), . . . , on voit sans peine 
que la fonction P(x) est négative dans des intervalles (-3, -2), (-5, -4), 
(-7, -6), . . . , c'est-à-dire dans tous les intervalles dont la limite inférieure 
est impaire. On peut donc passer tout de suite à l'intervalle (-4, -3) ; on 
trouve alors 

x + (x - l)(x - 2) + (x - l)(x - 2)(x - 3) + eP(x) 



eP(x-4) = 



(x - l)(x - 2){x - 3)(x - 4) 



Le dénominateur est positif, x étant compris entre et 1 ; quant an nu- 
mérateur, il est égal à 

x + (x-l)(x-2) 2 + eP(x), 
et il a le même signe que l'expression 

x 2 - x(l - x)(x - 2) 2 + exP(x); 

l 
mais le produit x(l - x) a pour valeur maximum -, puisque la somme de 

ses facteurs est égale à l'unité ; par conséquent, en remplaçant x 2 par 0, 
et (x - 2) 2 par 4, l'expression précédente est supérieure à 

ex P(x) - 1, 
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nombre évidemment positif. Ainsi la fonction P(x) est positive dans l'in- 
tervalle (-4, -3). 

Il reste à étudier le signe de P(x) dans des intervalles tels que (-6, -5), 
(-8, -7), . . . dont la limite inférieure est paire. Or, dans le domaine d'un 
pôle x - -m, on a 



P(x) = 



(-1) 



m 



1 • 2- ■■m x + m 



+ F(x + m), 



en désignant par F(x + m) une série entière enx + tn; par suite, si o est un 
nombre positif arbitrairement petit, pour x - -2n + o et x - —2n + 1 - o, 
on trouve 



P(-2n + o) 
P(-2n + 1 - o) = 



1 • 2 • • • 2n a 
1 



+ F(o), 
1 



l-2---(2n-l) a 



+ F(-o), 



de sorte que P(x) est positif à droite d'un pôle pair et à gauche d'un pôle 
impair. D'ailleurs, dans les intervalles considérés, P(x) peut admettre des 
valeurs négatives. En effet, par exemple, 



eP 



i-'-è) 



eP 



257 



eP 



755 
~32 



5 7 9 11 
2 ' 2 ' 2 ' ~2 



1 3 5 7 9 11' 

2 ' 2 ' 2 ' 2 ' 2 ' Y 



et comme 



( 



eP 



1 ! 1 

1 + TT + « ?r + 

3 3 5 



2 2 2 
égalité dont le second membre est inférieur à 2e, il s'ensuit que 



■(4) 



est négatif. On verrait de même que 



'■-y)' p (-îi 



le sont aussi. La fonction P(x) étant continue à l'intérieur de tout inter- 
valle limité par deux pôles consécutifs, on conclut de ce qui précède que 
l'équation P(x) - admet au moins une racine dans chacun des inter- 
valles f-y, -5j, f-6, - — \ f-y, -7), (-8, -y) , . . . tandis qu'elle n'( 



en a 
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point dans les autres intervalles. Si x est une valeur de la variable ap- 
partenant à un intervalle de limite inférieure paire et telle que P(x) soit 
négatif, comme P(x + 1) est aussi négatif, la relation 

P ( X ) ~ — 



entraîne l'inégalité 



d'où 



ex x 



P(x) - — > 0, 
ex 



TOI < -fr; 

ex 



la valeur absolue de P(x) correspondant à la valeur considérée de x est 
donc toujours assez petite. 

Ces résultats ont été obtenus par Bourguet ( x ) qui, en outre, a cherché 
à établir ( 2 ) que l'équation P(x) - ne peut avoir plus de quatre racines 
imaginaires, mais sa démonstration n'est pas rigoureuse. On n'a pas 
réussi jusqu'à présent à voir si l'équation P(x) - 0, dans les intervalles où 
elle possède une racine, n'en admet pas d'autres. L'étude de l'équation 
Q(x) - est encore moins avancée, car on ne sait même pas si elle possède 
des racines. Cependant Lindhagen ( 3 ) a prouvé qu'elle n'en a point si la 
partie réelle de x - 1 est négative. 



(!) Acta mathematica, t. II, 1883, p. 296-298. 

( 2 ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. XCVI, 1883, 
p. 1487-1490. 

( ) Studier ôfver gamma- funktionen, p. 34-37. 
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III. 

PROPRIÉTÉS DE LA FONCTION GAMMA. 



Relation des compléments. — Si l'on considère les deux expressions 



U(x) \ 1/ V n 



n + 1 -x _ r / x\ I x 

M1 + ----1 + 



n(l - x) n \ 1 / V n 

en les multipliant, on trouve 

1 n + l-x ( x 2 \ ( x 2 \ 



n(x)n(i-x) n \ i2j \ n 2 

mais 

sinra = ra (i-y(i-|)(i-|)..., 

par suite, à la limite, 

r(x)r(l-x) = -^—. 

sm 7TX 

Cette formule a reçu le nom de relation des compléments ; elle a été donnée 
par Euler dans ses Institutiones Calculi integralisi 1 ). Il est à peine besoin 
de faire remarquer que l'emploi de la formule de Weierstrass y conduit 
d'une façon immédiate. 

La relation des compléments, quand x est réel, permet de limiter 

l'intervalle où il est nécessaire de calculer T(x) à l'intervalle (0,-1. Dans 



( x ) 3 e éd., 1. 1, § 351, p. 218-221, et t. IV, p. 104-105. 
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le cas où la variable est imaginaire, le calcul de T(x), pour une valeur 

quelconque de x, pourra être ramené à celui de T(x) pour une valeur de x 

l 
dont la partie réelle soit comprise entre et -. 

I 
Si l'on fait x = - dans la relation des compléments, on obtient 



r(i)-* 



en prenant la détermination positive du radical, puisque la fonction T(x) 
est positive en même temps que x. Le résultat précédent, qui a une grande 
importance, est dû à Euler ( x ). 

Formule de Legendre. — Si dans le produit 

,-r, s Y 1 -2---W 

n(x) - n x 



x(x + 1) • •• (x + ri)' 

1 
on remplace x par x + -, on trouve 

m 1\ 2n x _l 2-4---2M 

fi \x+- n 2 



2/ 2x + 2n + l (2x + l)(2x + 3)---(2x + 2n-l)' 

d'où 

rr^rr^M 2n n 2x-l (2-4---2n) 2 

n(x)n x+- 2n 2- 



2/ 2x + 2n + l 2x(2x + l)(2x + 2)---(2x + 2n)' 

égalité dont le second membre peut se mettre sous la forme 

2 -2x 2n 2x l-2---2n (2-4---2n) 2 J_ 

2x + 2n + 1 ^ ; 2x(2x + 1) • • • (2x + 2ri) 1 • 2 • 3 • • • 2n yfc' 

Quand n augmente indéfiniment, le rapport 

2n 
2x + 2n + l 

tend vers l'unité, et l'expression 

,„ n2 x 1 • 2 • • • 2n 



2x(2x + 1) ■ ■ ■ (2x + 2n)' 



( ) Commentarii Academix Scientiarum imperialis petropolitanse, t. V, 1730-1731, p. 44. 
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que l'on obtiendrait en substituant 2x à x, et 2n à n dans Tl(x), devient 
égale à T(2x) ; enfin, le dernier facteur 

(2-4-6---2n) 2 2 
l-2-3---2n ^ 

pouvant être mis sous la forme 



, 2n + 1 / 2 • 4 • 6 • • • 2n \ 2 1 
2\ 2- 



2w \l-3-5---2n-l/ 2w + 1, 

on voit que sa limite, pour n - oo, est 2 V?? d'après la formule de Wallis ; 
de là résulte 

r(x)r ( X+ 2)^ r(2x) - 

Cette formule, trouvée par Legendre( 1 ), donne le moyen de réduire 

l'intervalle où il est nécessaire de calculer T(x) à l'intervalle 1 0, - J. Si x est 

imaginaire, la réduction porte naturellement sur la partie réelle de cette 
variable. 

La formule de Legendre n'est qu'un cas particulier d'une formule 
plus générale que nous allons donner. 

Formule de Gauss. — On a 

^/ s n r-1 1-2-n 

n(z) n x J 



n + x x(x + l)---(x + n - 1)' 
de même 

U[x+—)= : n x+ m i - — - 



\x+ — \[x+ — 

\ m ) \ m 



x+-+n \ x + -\\x+ - + 1 ... \x+ - + n-l 

m \ m I \ m j \ m 



m- 1\ n m-l_-| 

1I|X+ = : n X+ m i 

m I m- 1 
m ' x + + n 

m 

1-2---H 

x - 



/ m-l\( m-\ ^\ ( n-1 * 

\x+ x+ + 1 ••• \x+ + n - 1 

V m \ m \ m 



(- 1 ) Mémoires de la classe des Sciences mathématiques et physiques de l'Institut de France, 
année 1809, p. 485. Binet a donné de cette formule une démonstration directe différente 
de celle qui précède (Journal de l'École royale polytechnique, 27 e cahier, t. XVI, 1839, p. 209- 
210). 
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si l'on multiplie toutes ces relations, P désignant le produit des premiers 
membres, on obtient 



P = A- 
en posant 

A = 



mn mx- 
-n 



m+1 



(l-2---n) 



m 



mx + mn 



mx{mx + !)••• {mx + mn - 1) 



-m 



mn 



{mn) 



m-\ 



(mx + mn + l)(mx + mn + 2) • • • {mx + mn + m - 1)' 

fraction dont la limite est l'unité pour n = oo. 

D'autre part, si l'on remplace dans U(x) le nombre n par mn, et x 
par mx, on trouve une nouvelle expression 



Q = {mn) mx J 



1 • 2 • • • mn 



mx + mn ' mx{mx + 1) • • • (mx + mn - 1) 

dont la limite est T(mx) pour n - oo ; par suite 

.(l-2---n) m m mM+1 



P - AQm 



-mx' 



Soit 



d'où 



<p(ri) 



1 • 2 • • • mn ULiL 
n 2 



l-2---n 

ï ' 

n 2 



A — - — -m 2; 



Q (p(mn) 
quand n croît indéfiniment, le rapport 



<p m (n) 
(p(mn) 



a une limite A m qui est 



évidemment égale à celle de la fraction 

(p 2m (n) (p m (2n) 
(p 2 (mn) ' (p(2mn)' 



ou 



cp 2 (n) 
(p(2n) 



-\m v 



cp 2 (mn) 
cp(2mn) 



quotient dont les deux termes ont respectivement pour limites À? 2 et Aj ; 
or, 

(p 2 (n) __ p2(2-A-6---2n) 2 

cp(2n) ~ V n l-2-3---2n ' 



d'où 
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(p z (n) 2n + l( 2-4-6---2n 



(p{2n) \ 2n 11 • 3 -5- • -2n - 1/ 2n-l' 

par conséquent, d'après la formule de Wallis, À2 est égal à V2tt; on en 

conclut 

m—\ 
A m = (2rc) 2 , 

donc 

r(x)r xh — ■•■ru (271 ) 2 m 2 r(mx). 

Cette formule est due à Gauss ( x ) ; on doit y prendre les déterminations 
positives des radicaux, puisque le premier membre est positif pour toute 

valeur positive de x. 

l 

Si l'on pose x = — dans la formule de Gauss, elle devient 

1 m 

m—\ 



\tni \m) \ m / ^fm 

relation très élégante qui a été donnée par Euler( 2 ); on peut l'établir 
directement et en déduire la formule de Gauss par un procédé inverse 
de celui que nous avons suivi. 

Si maintenant on fait m - 2, on retrouve la formule de Legendre. Il 
est d'ailleurs facile, ainsi que l'a montré Binet ( 3 ), de tirer la formule de 
Gauss de la formule de Legendre. En effet, on vient de voir que le rapport 



U 



r(x)r(* + Il--r(* + — 

\ m) \ m 



m mx T(mx) 

est indépendant de x ; pour déterminer la constante U, il suffit de multi- 
plier l'expression précédente par celle que l'on obtient en y remplaçant 

1 

x par x H ; on trouve ainsi 

r 2m 

mT l x+ ±_)... T l x+ ^) T l x+ J_)... T l x+ 2 Jlzl 

2 \ m I \ m I \ 2m / \ 2m 

1 
m mx 2 Y{mx)l[mx + - 



( 1 ) Cari Friedrich Gauss Werke, t. III, p. 149-150. 

( 2 ) Opéra postuma Leonhardi Euleri mathematica et physica, t. I, p. 408-438. Voir aussi 
Bulletin des Sciences mathématiques, 2 e série, t. IV, 1880, p. 207-256. 

( 3 ) Journal de l'École royale polytechnique, 27 e cahier, t. XVI, 1839, p. 208-212. Parmi les 
nombreuses démonstrations de la formule de Gauss, celle de Sonine mérite aussi d'être 
citée (Bulletin de la Société mathématique de France, t. IX, 1880-1881, p. 162-166). 
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Les arguments des 2m facteurs du numérateur forment la suite croissante 



x, x + 
x + 



1 

2m' 
m 

2m' 



x + 



2m 
m + 1 



x + 



2m 



x + 



2m 



x + 



x + 

2m -1 
2m 



m 



2m 



de sorte que l'on peut grouper deux à deux ces facteurs de la façon 
suivante : 



r ^ r hl)}H x+ èù r hèA) 



x 



m — 1\ / m — 1 1 
2m I \ 2m 2 



en appliquant à chaque groupe la formule de Legendre, on a 



(2 yfnf 



c'est-à-dire 



r r ( 2x)r(2x + l)---r(2x 

{lyfnf 1 



+ 



m-\ 



m 



Mm 



-2mx 



^ ) 2mx+ m 2 



- T(2mx); 



d'autre part, en appliquant la formule de Legendre au dénominateur 

de U 2 , il devient 

1 



/// 



-2mx-k 2 V^ 
^2mx 



T(2mx); 



on en conclut que 



m-\ 



U = m2(2n) 2 , 

et la formule de Gauss est démontrée. 

On peut, au moyen de la formule de Gauss, en y donnant successive- 
ment à m les valeurs 3, 5, 7, 11, . . . restreindre de plus en plus l'intervalle 
où le calcul direct de T(x) est indispensable. On est ainsi conduit à se 
demander quel est l'intervalle minimum auquel on puisse parvenir ; Le- 
gendre ( x ) et, plus tard, Hoppe ( 2 ), se sont occupés de cette question, mais 
sans la traiter dans toute sa généralité. 

Formule de Mellin. — Soient 



«1/ «2/ 



a 



m 



(^Exercices de Calcul intégral sur divers ordres de transcendantes et sur les quadratures, 
1. 1, p. 283-285, et t. II, p. 26-32. 

( ) Journal fur die reine und angewandte Mathematik, t. XL, 1850, p. 152-154. 
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les m racines de l'équation binôme 

a m - 1 = 0; 



on a 



de même 



z - x - (z - a\x){z - a^x) • • • (z - a m x), 



;» 



X\ 



m 



nl \nl \ n nl\ n n 



1 + 



z x\ 

0Cm ~ I 
n ni 



enfin 

Si donc on pose 



Ot\ + «2 + " ' + a m - 0. 



P - T(z - a\x) T(z - aix) ■ ■ ■ T(z - a m x), 

en se reportant à la formule de Weierstrass relative à Fc(x) (p. 8), on 

trouve 

\ t— w y/ -r\tn /y\m-\/ M \mz 



\-r-*m+ïï-&YMïr- 



c'est-à-dire 



P= Z 



m 



i-|ï 



m 



n=l 



n=co 

n 

M = l 



1 + - 



z\ m I n smz 



ni \n + 1 



1- 



n + z 



m 



ou bien 



^n(-mrn[H^r 



par suite 



r m (z) 



n=oo 



T(z - a\x) T(z - a-ix) ■ ■ ■ T(z - a m x) 



n=0 



J- J- L \n + z 



m 



On doit cette formule à Mellin( 1 ). Mais Liouville( 2 ) a signalé le pre- 
mier l'intérêt que présentent des combinaisons de la forme 

T(aix) T((X2x) ■ ■ ■ T(a m x), 



C)Acta mathematica, t. III, 1883, p. 102-104. 

( 2 ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. XXXV, 1852, 
p. 321. 
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et il a même fait connaître la relation 



w 1+ a( 1+ a(i + e 



T(a 1 x)T(a 2 x)T(a 3 x) \ 1^)\ 2 3 )\ 3 3 

qui n'est qu'un cas particulier de la formule de Mellin. En effet, pour 
z—1, cette dernière formule se réduit à 



= n 



i- - 



r(i - a\x) r(i - «2 x ) ■ ■ ■ r(i - a m x) 

n=l 

et, si l'on fait m - 3, on obtient la relation indiquée par Liouville, tandis 
que, pour m - 2, on retrouve la relation des compléments. 

La formule de Mellin est susceptible d'être généralisée, comme Mellin 
l'a montré lui-même ( x ). 

Relation entre la fonction gamma et la série hypergéométrique. 

— La relation bien connue ( 2 ) 

y[(7 -l)-(2y-a-p- l)x] F + (y- a){y - fi)xF{y + 1) 

-y(y-l)(l-x)F(y-l) = 0, 

qui existe entre la fonction F et les deux fonctions contiguës F (y + 1) et 
F (y - 1), se réduit pour x - 1 à l'égalité 

F(a,jS,y,l) (y-«)(y-0) 



de même 



d'où 



F(a,p,y + 1,1) y(y-a-p)' 

F{a,§,y + 1,1) _ {y-a + l){y-§ + l) 
F(a,p,y + 2,1) ~ (y + l)(y -a -jS + 1)' 

F(a, j6, y + n - 1, 1) (y - a + n- l)(y - j8 + n - 1) 
F(a, f>,y + n,l) (y + n - l)(y - a - f> + n - 1) ' 

F(a, jS, y, 1) (y - a)(y - a + 1) ■ ■ ■ (y - a + n - 1) 



F{a,fi,y + n,l) y(y + l)---(y + n-l) 

(y-j8)(y-/3 + l).--(y-/3 + n-l) 



x 



(y - a - jS)(y -a-f> + l)---{y-a-f> + n-l)' 



( x ) Acta mathematica, t. III, 1883, p. 322-324. 

( 2 ) Cari Friedrich Gauss Werke, t. III, p. 130. Voir Maurice Godefroy, Théorie élémentaire 
des séries, p. 87-88. 
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c'est-à-dire, après avoir introduit au numérateur et au dénominateur les 
facteurs convenables, 

F(a, p, y, 1) {y + n)(y - a - p + n) U{y) U{y -a-§) 



F(a, f>, y + n, 1) (y - a + n)(y - p + n) Yl(y - a) Yl(y - p) 

Lorsque n augmente indéfiniment, F(a,p,y + n, 1) ayant pour limite 
l'unité, on trouve 

T(y)T(y - a - fi /U 

F(a,p,y,l) ^ 7- — — C . 

T(y - a) T(y - ]8) 

Les paramètres a, j6, y supposés réels doivent alors vérifier l'inégalité 

a + p -y < 0. 

La relation précédente peut servir, et c'est la méthode qu'a suivie 
Gauss, à établir les principales propriétés de la fonction gamma. C'est 
ainsi qu'il en déduit, par exemple, la relation des compléments de la 
manière suivante : 

On a 



,'n n 1 



r 2/l 



2 r r i r /i _«) r /i + « 



or, 



in n 1 \_ n 2 n 2 (n 2 -2 2 ) n 2 (n 2 - 2 2 )(n 2 - 4 2 ) 

\2'~2'2' / ~ 1~ 2 + 1-2-3-4 1-2-3-4-5-6 + '"' 

série que l'on obtient en faisant x = - dans le développement de cos nx 
en fonction de sin x : 

n 2 . 2 n 2 (n 2 -2 2 ) . 4 n 2 (n 2 - 2 2 )(n 2 - 4 2 ) 



2 n. in — t. i . a n in — t. )\n — •* ) . f, 

sin x + — — - — - — — sin x — - — - — - — - — - — sin x + 



1-2 1-2-3-4 1-2 -3 -4-5 -6 

par suite 



Z A! \A Al cos 

2 

ou bien, en remplaçant n par 2x 



2 / \2 / cos7ix 



( 1 ) Cari Friedrich Gauss Werke, t. III, p. 147. La généralisation de cette formule a fait 
l'objet d'un Mémoire de Mellin (Acta Societatis Scientiarum fennicx, t. XXIII, 1897, n° 7). 
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l 
ou encore, par la substitution de x à - - x, 

T(x)T(l-x)= -^—. 
sm7ix 



Ainsi la théorie de la fonction gamma se déduit sans difficulté de celle 
de la série hypergéométrique. D'ailleurs inversement, comme l'a montré 
Thomae ( x ), la théorie de la série hypergéométrique peut être exposée, 
d'une façon tout à fait élémentaire, en prenant pour point de départ les 
premières propriétés de T(x). 



(*) Zeitschrift fur Mathematik und Physik, t. XXVI, 1881, p. 314-333, et t. XXVII, 1882, 
p. 41-56. 
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IV. 

FONCTION DE BINET. 



Soit 



L'inégalité 



K„4 + ^±)log(l + ^-)-l. 



\x + n\ > 1 
est vérifiée, quel que soit x, à partir d'une certaine valeur de n ; alors 

lo (i , 1 \ = 1 1 1 , 1 : 

° g \ + x + nl x + n 2 ( x + n) 2 + 3 (x + n) 3 '"' 

et l'on peut poser 



Mm = 



(x + n) z 

en désignant par A n un coefficient qui reste fini quand n augmente au 
delà de toute limite ; par suite, si r est le module de x, et À une constante 
positive supérieure, quel que soit n, au module de A n , on a, à partir d'une 
valeur suffisamment grande de n, 

, , A 



< ~ \2' 



(n - r) z 

la série de terme général u n est donc absolument et uniformément conver- 
gente pour toute valeur de x autre que zéro ou qu'un entier négatif. La 
somme de cette série se représente par cd(x) ( x ), de sorte que 

M=0 



() C'est Cauchy qui le premier a fait usage de cette notation (Œuvres complètes, 
l re série, t. VII, p. 281). 
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La fonction cd(x), introduite par Plana ( x ), est généralement appelée fonc- 
tion de Binet. On doit, en effet, à Binet( 2 ) une étude très développée de 
cette transcendante. 

Si l'on considère les deux développements 



[ x+ l)H i+ D 


-i 


+ 


[ x+ î)M i+ , + i) 


-1 


[( X+ 2M 1 + * + l) 


-i 


+ 


[(x + -)log(l + i + 2 ) 


-1 



+ 



+ 



on obtient par soustraction 

cd(x) - <d(x + 1) = lx+ -jlog(l + -J-l; 

cette relation, signalée par Binet ( 3 ), joue, par rapport à la fonction cd(x), 
un rôle analogue à celui de la relation fonctionnelle dans l'étude de T(x). 
Lorsque n augmente indéfiniment, cd(x + n) tend vers zéro, comme il 
est facile de s'en convaincre en se reportant à la série qui définit cd(x). Ce 
résultat, joint à la relation fonctionnelle qui vient d'être établie, suffit à 
caractériser la fonction de Binet. En effet, de l'égalité 

<D(X) - (D(X + 1) = \X + - jlog (l + -) - 1/ 

on déduit successivement 

cd(x + 1) - o(x + 2) = [x + -) log (l + — ) - 1, 

cù(x + n) - cd(x + n + 1) - lx+ )l°g(l + ) ~ 1/ 



d'où 



(D 



Wf/ 2p + l\ ( 1 

(x) - a(x + n + 1) = 2^ k+— ^— llogll + 



X + p t 



-1 



et par suite, pour n - oo, 



<D(X) - Yj 



n=0 



X + 



2n + V 



W 1 *^;)- 1 



( ) Memorie délia reale Accademia délie Scienze di Torino, t. XXIV, 1820, p. 410. 
() Journal de l'École royale polytechnique, 27 e cahier, t. XVI, 1839, p. 220-269. 
( ) Journal de l'École royale polytechnique, 27 e cahier, t. XVI, 1839, p. 228. 
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Formules de Binet. — Le terme de rang n de la série o(x) peut se 
mettre sous la forme 

« n _i = -(x + n) log 1 1 + - log ( 1 ) - 1; 

n i & \ x + n/ 2 & \ x + n/ 

mais, pour toute valeur entière de n, 

I 1 \_ 1 11 11 11 

° S \ x + n)~ x + n 2 (x + n) 2 3 (x + n) 3 4 (x + n) 4 '"' 

en supposant qu'à partir de n - 1 l'inégalité 

|x + n\ > 1 

soit vérifiée, c'est-à-dire que le point x reste à l'extérieur d'une série de 
cercles de même rayon égal à l'unité décrits des points -1, -2, -3, . . . 
comme centres ; le n ieme terme de cd(x) est alors développable suivant la 
série absolument convergente 

1 1 2 1 3 1 

""-! " 3^4 (x + n) 2 + 4~ 6 (x + n) 3 + 5^8 (x + n) 4 + "' 

de plus, si l'on pose 

1 12 13 1 

3 • 4 |x + n\ 2 4 • 6 \% + n\ 3 5 • 8 \% + n\ 4 

on en déduit 

1 1 

12 |x + n| (|x + n\ - 1) 

inégalité dont le second membre est le terme général d'une série positive 
convergente ; par suite, d'après le théorème des séries de séries, 



n=co n=co 

( \ - — — V —— V 

a[X) " 3 . 4 Zj < x + n \2 + 4 . 6 Li (x + n) 3 



n=co 

3 V 1 ! i i -, 

+ -^^ / 7 + -"/ x + n > 1. 

5 • 8 ^-j (x + n) 4 



Cette formule remarquable a été établie pour la première fois par Binet ( 1 ). 



( 1 ) Journal de l'École royale polytechnique, 27 e cahier, t. XVI, 1839, p. 226. La démons- 
tration que nous donnons ici nous parait plus simple que toutes celles qui ont été 
proposées. Voir, par exemple, celle exposée par Bourguet dans sa Thèse, Annales scien- 
tifiques de l'École normale supérieure, 2 e série, t. X, 1881, p. 178-185. 
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Maintenant le ternie de rang n + 1 de cd(x) peut être mis sous la forme 

u n = (x + n)log(l + |+ -logll + )- 1; 

v ; 6 \ x + n) 2 6 \ x + n/ 

or, à partir de n = 0, on a 

/ 1 \ _ 1 11 11 11 

0g l + x + n) ~ x + n ~ 2 (x + n) 2 + 3 (x + nf ~ 4 (x + n) 4 + " ' ' 

en supposant que, pour toute valeur nulle ou entière de n, l'inégalité 

|x + n\ > 1 

soit satisfaite, c'est-à-dire que le point x reste à l'extérieur d'une série de 
cercles de même rayon égal à l'unité décrits des points 0, -1, -2, -3, . . . 
comme centres ; par conséquent 

11 2 1 3 1 

u n = - — - - — r + 



3 • 4 (x + n) 2 4 -6 (x + n) 3 5 • 8 (x + n) 4 

d'ailleurs, la série positive de terme général 

11 2 1 3 1 

v n = 7T- Â -, ^ + 1—7: ^ + 



3-4|x + n| 2 4-6|x + n| 3 5-8|x + n| 4 
est convergente, comme on vient de le voir ; donc 



= J_ y 1 _ J_ y 1 



3 • 4 -t_j (x + n) 2 4-6 ^— ' (x + n) 3 

n=co 

3 Y^ - 1 11-1 

+ ^^T / 7 ---- / x + n > 1, 

5 • 8 Lu (x + n) 4 

formule due également à Binet ( x ). Ce dernier développement diffère du 
précédent par l'alternance des signes; en outre, dans les sommes, la 
valeur initiale de n est zéro au lieu d'être l'unité. 

L'étude des séries de Binet, leur transformation en séries se prêtant 
mieux au calcul numérique, et enfin la détermination du reste, ont donné 



( 1 ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. IX, 1839, p. 158. 
On a souvent attribué à Féaux la seconde des formules de Binet. 
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lieu à de nombreux travaux dont les principaux sont ceux de Genocchi ( 1 ), 
Gilbert ( 2 ), de Tilly ( 3 ), et Bourguet ( 4 ). 

Formule de Gudermann. — Si au développement 

logll(x) = xlogn-logx- logfl + -J- logfl + -J logfl + -J 

on ajoute les deux termes 

xlog(l + ï) + log(l + ï) 

qui tendent vers zéro pour n - oo, le nouveau développement obtenu 

u n = xlog(n + x)-logx-logfl + -J logfl + — — -J 

a également pour limite log T(x) quand n devient infini. De cette dernière 
relation on déduit 

U n+1 -U n = xlog(l + J- logfl + -J, 

\ x + ni \ ni 

différence dont le second membre peut se mettre sous la forme 



.( x+ ^) log ( 1 + xT^) 



-î 



2n + l 



log(l + i)-l 



2n + l 



log 1 + 



n + 1 



2 

2n - 1 , / x 

— ^ lo g 1 + - 



or 



M = l 

et, pour x - 1, on trouve 



n=\ 



( ) Bulletins de l'Académie royale des Sciences, des Lettres et des Beaux-Arts de Belgique, 
t. XX, 1853, p. 392-397, et t. XXI, 1854, p. 84-95; 2 e série, t. XXXVI, 1873, p. 546-565, et 
t. XXXVII, 1874, p. 351-352. —Annali di Scienze matematiche efisiche, t. V, 1854, p. 150-160 ; 
2 e série, t. II, 1859, p. 380. 

( 2 ) Mémoires de l'Académie royale des Sciences, des Lettres et des Beaux-Arts de Belgique, 
in-4°, t. XLI, 1875-1876 (n° 3), p. 14-22. 

( 3 ) Bulletins de l'Académie royale des Sciences, des Lettres et des Beaux-Arts de Belgique, 
2 e série, t. XXXV, 1873, p. 30-40, et t. XXXVIII, 1874, p. 67-70. 

( ) Annales scientifiques de l'École normale supérieure, 2 e série, t. X, 1881, p. 178-191. 
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d'autre part, comme on le reconnaît facilement, 

Ï'Z — CO 

1, „ > V [ 2 " + l, /, 1 \ 2«-l, r x\ 
x- -log(l +,) = ^ [— log(l + —)- — log(l + ") 

M = l 

enfin 

M=oo 

log T(x) - X log(l + x) - log X + )| X log 1 H I - log 1 + - I 

par suite, en appelant À la différence 1 - cd(x), on a 

logT(x) - À + Ix- -Jlogx-x + <d(x). 

Pour déterminer la constante À, il suffit de donner successivement à x les 
valeurs n et 2n ; on obtient de cette manière 

2 log r(n) = 2A + (2n - 1) log n - 2n + 2<D(n), 

logT(2n) = A + (2n - - Jlog2n - 2n + o(2n), 

d'où, par soustraction, 

A - log A M + cD(2n) - 2cD(n), 

en posant 

1 on-1 U-2---(n-l)] 2 
A "^ 22 2 l-2---(2n-l) ' 

expression à laquelle on peut donner la forme suivante 

Ï2 (2 • 4 • 6 • • • 2n) 2 
An " Vn l-2-3-..2n ' 

ou bien 



, 2n + 1 / 2 ■ 4 • 6 ■ ■ ■ 2n \ 2 1 

A„ = a 4- 



2n \l-3-5---2n-l/ 2n + l' 

quand n augmente indéfiniment, A n a donc pour limite \2n d'après la 
formule de Wallis, et, comme cD(2n) et cD(n) s'annulent, on en conclut que 
la constante A est égale à log V2h. Ainsi 

logT(x) - log V27T + (x- -Jlogx-x + cD(x); 
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c'est la formule de Gudermann i 1 ). La démonstration précédente, qui re- 
vient au fond à transformer la définition même de T(x), est due à Jen- 
sen ( 2 ). 

Application. — On peut, au moyen de la formule de Gudermann, 
trouver une nouvelle expression du module de T(x) pour une valeur 
imaginaire de x. 

Soit, en effet, 

x - a + fti. 

On a 

? ? - P 

log(a + fti) - log(a + jS )2 + i arctang -, 

d'où 

1 1/ 1\ -P n P •/ 1\ P 

X X ~\ = [ a 2 + ^ ] 2( a -2J +l 2fP aTctan ëâ +l ( a -2J aTctan ëâ / 

par suite 

\ X X - 1 2\ = [a 2 + p 2 ]i( a ~V ^arctang §. 



d'autre part, 



|e- x | - e" a ; 



enfin, le module de cD(a + j3z) tend vers zéro quand a ou j3 croissent 
indéfiniment, car il vérifie l'inégalité 



n=oo 



|û>(a + /3i)| < ^2 £ 



n=0 J(â + n) 2 + |S 2 (J(â + n) 2 + jS 2 - l) 
on peut donc poser 

|r(a + /3i)| - V2^[a 2 + i S 2 ]2r y e -"-^ arctan §!(l + £ ), 

en désignant par e un nombre qui s'annule pour a - oo ou f> - oo. Cette 
formule est due à Mathias Lerch ( 3 ). 



() Journal fur die reine und angewandte Mathematik, t. XXIX, 1845, p. 210. Gudermann, 
dans son Mémoire, semble craindre qu'il n'y ait contradiction entre sa formule et la 
formule de Stirling. Cette appréhension est rien moins que justifiée, comme on va le 



voir. 

a 



( z ) Nyt Tidsskrift for Mathematik, t. II B, 1891, p. 40-41. 
( 3 ) Théorie funcke gamma (Véstnîk ceské Akademie cîsafe Frantiska Josefa pro vëdy, sloves- 
nost a umêni v Praze, t. II, 1893). 
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Formule de Stirling. — Soit x une variable réelle et positive. Le 
n ieme j- erme de la série cd(x) est développable suivant la série positive 
convergente 

11 2 1 3 1 

u n-l ~ JTg ( X + M )2 + F6 (x + n) 3 + 5^8 (x + n) 4 



et vérifie l'inégalité 



U M-1 < T7T 



12 (x + n - l)(x + n)' 



d'où l'on déduit 



1/1 1 \ 
M0< Î2lx-^Tj' 



1 /l 

12 \x 

1/1 1 



«1 < — 



12\x + l x + 2/' 



U M-1 < T7T 



12 \x + n - 1 x + n 



on peut donc poser 

o(x) = 



12*' 

en désignant par 6 un nombre positif compris entre et 1 . Si l'on remplace 
cd(x) par cette expression dans la formule de Gudermann 

logT(x) - log V27Ï+ (x- -Jlogx- x + ©(*), 

elle devient 

T(x) - V27ix x ~2e~ x+ Ï2^; 

c'est la formule de Stirling i 1 ). Lorsque la variable est égale à un entier m, 
cette formule se réduit à 

1 R 

l-2-3---m= y/2nm m+ 2e~ m+ T2E, 

résultat qui a une grande importance dans le calcul des probabilités. 

Au lieu de déduire la formule de Stirling de la formule de Guder- 
mann, on peut aussi l'établir directement, et d'une façon très simple, au 



(^Methodus differentialis : sive Tractatus de summatione et interpolatione serierum infini- 
tarum, p. 135-139. 



44 



moyen de la formule de Wallis. Ce dernier procédé a fait l'objet de bien 
des démonstrations, dont notamment celles de Glaisher( 1 ), Cayley( 2 ), 
Mansion ( 3 ), Cesàro ( 4 ), Rouché ( 5 ). 

Applications. — I. La formule de Stirling permet de calculer deux 
limites du nombre de Bernoulli B n . 

En effet, 

_ l-2-2n c 
Bn ~ 1 2n-l n 2n Sln ' 

par suite, 

Bn __ l-2---2n S 2n 

~Ë[ ~ 2 2M-l 7I 2n-2 ~~S^' 

mais S<in décroissant à mesure que n augmente, il en résulte 

B n 1 • 2 • • • 2n 

< 



B 1 2 2n - 1 n 2n - 2 ' 

d'où, comme B^ est égal à -, 

^ 1 T(2n + 1) 
" 12 (27i) 2 "- 2 ' 

d'autre part, si, dans l'expression de B n , on remplace Si n par l'unité, on 
en déduit 

r(2n + 1) 

Or, d'après la formule de Stirling, T(2n + 1) vérifie les inégalités 

T(2n + 1) > V2^(2n) 2M+ 2 e - 2M 7 

2 1 

T(2n + 1) < V27z(2n) 2M+ 2e" 2M+ 24^, 



(- 1 ) The quarterly Journal ofpure and applied Mathematics, Vol. XV, 1878, p. 57-63. 

( ) The quarterly Journal ofpure and applied Mathematics, Vol. XV, 1878, p. 63-64. 

( ) Nouvelle Correspondance mathématique, t. V, 1879, p. 44-51. 

( ) Nouvelle Correspondance mathématique, t. VI, 1880, p. 354-357. 

( ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. CX, 1890, 
p. 513-515. La démonstration de Rouché, comme toutes les précédentes, est relative au 
cas où x est un entier positif. Gomes Teixeira a étendu cette démonstration au cas où 
x est un nombre positif quelconque (Nouvelles Annales de Mathématiques, 3 e série, t. X, 
1891, p. 312-315). 
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on a donc 

1 (2nf n+ i -2n+^r- 

Bn< Ï2 S 6 ' 

(2n) 2n ~2 

(2n) 2n+ 2 , 
B n >2 y —t -e- 2n , 

(2n) 2n ~2 

La première de ces limites, entre lesquelles B n reste compris, a été 
indiquée par Cauchy ( 1 ). 

IL La formule de Stirling est d'un usage commode pour l'évaluation 
de certaines limites. 

Soit, par exemple, 

_ (1 • 2 • • • n) m m mn+1 m 

n ~ 1 • 2 • • • mn muL ' 

n 2 



on a 



m , m 



par suite, 



(1 • 2 • • • n) m = (2n)2n mn+ 2 e - mn+m0 ( n \ 

1 1 

1 • z • • • mn — (In) ^ {mn) ^ e v ', 



1 Til — 1 

U n = m2(2n)^e m ^ n) -^ mn ^; 



quand n croît indéfiniment, m restant fixe, a)(n) et oimri) tendent vers 
zéro, et U n a pour limite 



m—\ 



U = m2(2n) 2 . 

C'est ainsi que Cauchy ( 2 ) a déduit la formule de Gauss de la formule de 
Stirling. 



( ) Exercices d'Analyse et de Physique mathématique, t. II, p. 396-397. 

( ) Exercices de Mathématiques (Œuvres complètes, 2 e série, t. VII, p. 121-123). 
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V. 



FONCTIONS O(x) ET W(x). 



On a, pour toute valeur de x non égale à un entier négatif, 



logT(x + 1) - -px + 



f-log(l + f)] + [f-log(l + f 



+ • 



Or, soit r un entier positif arbitrairement grand ; à partir de n = r + 1, 
l'inégalité 



n(x + n) 



< 



n(n - r) 



est vérifiée si la variable x reste à l'intérieur du cercle de rayon r dé- 
crit de l'origine comme centre; la série logT(x + 1) est donc dérivable, 
puisque la série formée par les dérivées de ses termes est uniformément 
convergente. Ainsi 

d\ogT(x + l) i 1 \ /l 1 \ /l 1 

- = -p+ ' 



1- 



+ -- 



+ -- 



+ ... 



dx r ' \ x + U ' \2 x + 2j ' \3 x + 3, 

Ce développement était connu d'Euler ( x ). Il offre la particularité re- 
marquable que la série des fractions simples qui y figurent devient ab- 
solument convergente par l'addition d'une constante à chacune de ces 
fractions. Ce fait, d'après Weierstrass, a été la première indication qui ait 
mis sur la voie du théorème de Mittag-Leffler ( 2 ). 

On obtient de même, pour toute valeur de x non égale à un entier 
négatif, 

d 2 logT(x + l) 



dx 2 



+ 



+ 



(x + iy (x + iy (x + 3y 



+ 



( ) Institutiones Calculi differentialis, pars posterior, § 384, p. 801. 

( ) Voir une lettre adressée par Hermite à Mittag-Leffler et reproduite dans le Journal 
fur die reine und angewandte Mathematik, t. XCII, 1882, p. 145-155. 
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développement dû à Gauss ( x ) ; en dérivant encore une fois, on trouverait 
une nouvelle série, et ainsi de suite. En général 



(-1)" d n \ogT(x + l) 1 1 1 

+ Z T7-T + Z T7-T + 



1-2- ••(«-!) dx n (x + l) n (x + 2) n (x + 3) M 

Nous poserons 

d log T(x) d 2 loe T(x) 

d'après ce qui précède, pour toute valeur de x autre que zéro ou qu'un 
entier négatif, ces fonctions, sommes de séries absolument et uniformé- 
ment convergentes, sont holomorphes. 

Les fonctions O(x) et ^(x), que l'on a désignées par les notations 
les plus diverses, jouent un rôle important dans la théorie de la fonc- 
tion gamma. C'est ainsi que Hôlder ( 2 ) est parvenu à établir que T(x) 
n'était solution d'aucune équation différentielle algébrique, en prouvant 
au préalable que O(x) avait ce caractère. La fonction W(x) ne présente pas 
moins d'intérêt, car, en partant de son développement, on peut retrouver, 
comme l'a montré Hermite ( 3 ), toutes les propriétés de T(x). Inversement, 
il est facile de déduire toutes les propriétés de O(x) et de W(x) de celles 
de T(x) ; si l'on prend, en effet, les dérivées logarithmiques des relations 

T(x + 1) = xT(x), 
T(x + m) = x(x + 1) • • • (x + m - l)r(x), 

r(x)r(i-x) n 



smnx 

r(x)r xh — ---r xh \ = (2n) 2 m 2 Timx), 

\ m! \ m 1 



( ) Cari Friedrich Gauss Werke, t. III, p. 153. 

( ) Mathematische Annalen, t. XXVIII, 1887, p. 1-13. — La fonction gamma est, parmi 
les fonctions transcendantes élémentaires, la seule qui présente cette particularité. Au 
sujet du théorème de Hôlder, on pourra consulter un Mémoire de Barnes inséré dans 
The Messenger of Mathematics, t. XXIX, 1899, p. 122-128. 

( 3 ) Faculté des Sciences de Paris. — Cours de M. Hermite rédigé en 1882 par M. Andoyer, 
3 e éd., p. 113-118. 
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on obtient 



1 
0(x + 1) - O(x) = -, 

x 

1 1 
0(x + m) - O(x) - - + + • • • + 



x x + 1 x + m - 1' 

O(x) - 0(1 - X) - -71 COt 7ZX, 

1 \ . / m — V 



O(x) + Oxh — ) + ••• + O I x H I - mO(mx) - m log m, 

et, en dérivant encore une fois, 



W(x + l)-W(x) = -\, 



x L 



^(x + m)-^(x) = -^r 



x 2 (x + l) 2 (x + m - l) 2 

W(*) + W(l -x) = l -^— ) , 

V sin tcx I 

W{x) + W (x +—) + ■■•+ W (x + 1 H—1\ = m 2 w{mx). 
\ m) \ m ! 

Enfin, à toutes ces formules on peut joindre les deux suivantes 

,l-2---n 



W {n \x + l)-W in) (x) = (-l) n - 



x n+1 



p=m 



n+1' 



-[ (X + p - 1) 

que l'on trouve en dérivant n fois les deux premières relations relatives 
à O(x). 

Nous avons vu (p. 7 et 9) que l'on a formé de nombreuses fonctions 
semblables à T(x) ; on a été amené, par là même, à construire des fonctions 
analogues à O(x) ou à ^(x). Telle est, entre autres, la fonction 

<p(x + ï) = — [xO(x + l)], 

étudiée par BlasernaC 1 ) à propos de recherches optiques, et dont les 
propriétés rappellent celles de O(x). 



( ) Memorie délia reale Accademia dei Lincei, 5 e série, 1. 1, 1896, p. 499-557. 
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/l l l \ 
Limite de l'expression log n — I - + + • • • + 1 pour n = oo. 

y A. JL \ L A. r fi L / 

— La série 

peut se mettre sous la forme 

(1 1 1 1 \ /l 1 1 \ 

\1 2 3 n r I \x x + 1 x + n-U 

en désignant par R n le reste ; or, si l'on pose 

11 1 , 

P«^ ï + 2 + '" + n~ gM ' 

la série précédente devient 

/Il 1 \ 

O(x) - p n -p + logn- [- + - + ••• + ) + R n ; 

r r \x x + 1 x + n-ll 

on en conclut que O(x) est la limite, pour n - oo, de l'expression 

logn-(- + - + ••• + -), 

\x x + 1 x + n-ll 

résultat qui pourrait servir à définir O(x) i 1 ). En particulier, pour x = 1, 
on a 

P = -$(1). 

Développement de la différence 0(x + a) - O(x). — De la formule 

on déduit que la différence 

O(x) - 0(x - a) 
est égale à 

/ 1 1\ / 1 1 \ / 1 1 \ 

\x-a x) \x-a + l x + 1/ \x-a + 2 x + 2) 



( 1 ) Voir un Mémoire publié par Eduard Weyr dans Casopis pro pëstovânî mathematiky 
afysiky, t. XXI, 1892, p. 151-166. 
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développement qui est susceptible d'être transformé de deux façons 
différentes. 

Si l'on suppose d'abord l'inégalité 

\a\ < \x + n\ 

vérifiée à partir de n - 0, les termes peuvent être considérés comme les 
sommes de séries entières absolument convergentes telles que 



+ 



+ 



x-a + n x + n (x + n) z (x + n) à (x + np 



+ 



Or, soient r le module de x, s celui de a, et m l'entier immédiatement 
supérieur à r + s ; à partir de n = m, on a 



I 

p=l 



aP 



(x + nf +l 



< 



E 



sP 



T[{n- rf + 



1' 



ou, en désignant par U n le premier membre de cette inégalité, 

1 1 



U n < 



n- r - s n - r 



il résulte de là que la série de terme général U n est convergente, puisque, 
à partir de n - m, ses termes sont inférieurs à ceux de la série positive 
convergente 

n=oo 

0(m - r) - 0(m -r - s) - / ) ; 

*—'\tn-r-s + n m-r + nl 
n=0 

on peut donc, d'après le théorème des séries de séries, sommer par 
colonnes les développements 



x-a 




x 


1 




1 



+ 



+ 



+ 



+ 



x-a + 1 x + 1 (x + 1) 2 (x + 1) 3 (x + 1) 



+ 



+ • 



+ 



+ 



x-a + n x + n (x + n) z (x + n) (x + np 



+ 
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par suite 

n=oo n=oo n=co 



O(x) - 0(x -a) - a; - + a 2 ; + a 3 ; + 

L -'(x + n) 2 L - i n (x + nf L ~'(x + n) 4 

ou, en remplaçant a par -a, 

n=co n=oo n=oo 
®(x + a) - O(x) = a ) - a 2 } + a 3 ) - 

n=0 v ' n=0 v ' n=0 v ' 



c'est-à-dire 



2 

ci a 
0(x + a) = O(x) + -O'(x) + 0"(x) + • 



Cette série, qui a été considérée par Gauss ( x ), s'obtiendrait tout de suite 
par l'application du théorème de Cauchy relatif au développement d'une 
fonction en série entière ; mais il nous a paru intéressant de suivre une 
marche identique à celle adoptée pour la seconde transformation. On 
doit supposer que le module de a est inférieur au plus petit des modules 
des nombres x, x + 1, x + 2, . . . , x + n, . . . , ou, autrement dit, que le point x 
reste à l'extérieur d'une série de cercles de même rayon égal au module 
de a décrits des points 0, -1, -2, -3, . . . comme centres; cette condition 
concorde avec celle donnée par le théorème de Cauchy. 

Si maintenant ou suppose positive la partie réelle de la différence 
x - a, les termes de la série 



\X-a x) \x - a + 1 x + 1/ \x-a + 2 x + 2, 
sont développables en séries absolument convergentes (p. 15) telles que 
1 1 



x - a + n x + n 

a a(a + 1) 

(x + n)(x + n + 1) (x + n)(x + n + l)(x + n + 2) 

d'ailleurs, soient encore r le module de x, s celui de a, et m l'entier 
immédiatement supérieur kr + s; k partir de n = m l'inégalité 

p=co 

a(a + 1) • • • (a + p - 1) 



L 



(x + n){x + n + 1) • • • (x + n + p) 



s(s + 1) • • • (s + p - 1) 



ES^S ï- L) • • • 
(n - r)(n - r + 1 



(n - r)(n - r + 1) • • • (n - r + p) 



( 1 ) Cari Friedrich Gauss Werke, t. III, p. 153. 
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est satisfaite, d'où, en représentant par V n la première somme, 

1 1 



V n < 



n-r - s n-r 



on en conclut, comme précédemment, que l'on peut sommer par co- 
lonnes les séries 

1 1 a a(a + 1) 

+ — Z-, T7 + --- , 



x - a x x(x + 1) x(x + l)(x + 2) 

1 1 a a(a + 1) 



x - 


-fl + 1 


x + 1 


(x + l)(x + 2) (x + l)(x + 2)(x + 3) 




1 


1 

x + n 


a a(a + 1) 


x - 


- a + n 


(x + n)(x + n + 1) (x + n)(x + n + l)(x + n + 2) 



mais 
1 



+ ^— 7T 7—. -7T + ---/ 



p x(x + 1) • • • (x + p - 1) x(x + 1) • • • (x + p) (x + 1) • • • (x + p + 1) 

par conséquent 

a 1 a(a + ï) 1 a(a + l)(a + 2) 
O(x) - 0(x - a) = - + - —. -^ + - -; -^7 rr + 

x 2 x(x + 1) 3 x(x + l)(x + 2) 

ou bien, en changeant le signe de a, 

, ^, s « 1 «(« + 1) 1 flffl + l)(fl + 2) 

0(X + fl) - O(X) = ; + ; ri ~ 

v ' w x 2 x(x + 1) 3 x(x + l)(x + 2) 



la partie réelle de x + a étant positive. 

En particulier, si a est égal à un entier m, cette formule devient 

^, . _ . m 1 m(w - 1) lm(w-l)(w-2) 

0(x + m) - O(x) = ; + ^ 

v ; w x 2 x(x + 1) 3 x(x + l)(x + 2) 

(-l) m w(m - 1) • • • 1 

m x(x + l)---(x + m - 1)' 

d'où l'identité 

11 1 ml m(m - 1) 

x x + 1 x + m-1 x 2 x(x + 1) 

1 m(m - l)(m - 2) (-l) m m(m -!)•••! 



+ 



3 x(x + l)(x + 2) m x(x + 1) • • • (x + m - 1)' 
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qui, pour x = 1, se réduit à la suivante : 

11 \ m 1 m(m - 1) lw(w-l)(w-2) . „.,„ 1 
- + - + ••• + — = i i + i - i (-1) — 

12 m 1 2 1-2 3 1-2-3 v 7 m 

relation due à Jean Bernoulli. 

Si l'on égale les deux développements dans lesquels on a successive- 
ment transformé la série 

O(x) - 0(x - a), 

on trouve 

n=oo n=oo n=°o 1 

a / ô + a 2 / ô + fl 3 / r + • • ■ 

*-■ (x + n) 2 L ~' n (x + nf L ~' n (x + n) 4 

a 1 a(a + 1) 1 a(a + l)(a + 2) 
— 1 h — + • • • 

x 2 x(x + 1) 3 x(x + ï)(x + 2) 

identité très remarquable que Binet ( x ) a utilisée dans son Mémoire sur 
les intégrales eulériennes ; pour qu'elle soit valable, il faut que a ait une 
partie réelle moindre que celle de x et un module inférieur au plus petit 

des modules des nombres x, x + 1, x + 2, . . . , x + n, 

L'identité précédente donne, lorsqu'on identifie les coefficients de a 
dans les deux membres, 

.111 1 1-2 

y ' x 2 x(x + 1) 3 x(x + l)(x + 2) 

en admettant que la partie réelle de x soit, ou positive, ou négative mais 
d'une valeur absolue moindre que le plus petit des modules des nombres 
x, x + 1, x + 2, Cette série a été étudiée par Bauer ( 2 ). 

Détermination de O(x) pour une valeur rationnelle de x. — La 

fonction O(x) est susceptible d'être mise sous forme finie toutes les fois 
que x est un nombre rationnel. On l'établit de la manière suivante, au 
moyen d'une méthode extrêmement ingénieuse qui a été indiquée par 
Gauss( 3 ). 



(^Journal de l'École royale polytechnique, 27 e cahier, t. XVI, 1839, p. 231 et p. 255. 

( 2 ) Journal fur die reine und angewandte Mathematik, t. LVII, 1860, p. 256-272. 

( 3 ) Cari Friedrich Gauss Werke, t. III, p. 155-157. Un autre procédé, plus rapide, a été 
donné par Jensen (Nyt Tidsskrift for Mathematik, t. II B, 1891, p. 52-54). Voir aussi, à ce 
sujet, un Mémoire de Mathias Lerch (Rozpravy ceské Akademie cîsafe Frantiska Josefa pro 
vëdy, slovesnost a umëniv Praze, t. V, 1896, n° 14, p. 19-21). 
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Si l'on suppose d'abord le nombre rationnel x extérieur à l'intervalle 
(0, 1), en désignant par m l'entier positif immédiatement inférieur à x ou 
supérieur à -x suivant que x est positif ou négatif, les deux formules 

/ 1 1 1 \ 

O(x) - 0(x - m) + + + • • • + . 

\x-l x-2 x-mJ 

(1 1 1 \ 
O(x) = 0(x + m)-\- + - + ••• + , 

\x x + 1 x + m-ll 

montrent que le calcul de O(x), pour une valeur rationnelle quelconque 
de x, peut toujours être ramené au calcul de O(x) pour une valeur ration- 
nelle de x comprise entre et 1. 

Soit donc 

V 
x- -, 

q 

l'entier p étant moindre que l'entier q. On a vu (p. 49) que O(x) peut se 
mettre sous la forme 

/Il 1 \ 

O(x) = p n - p + log n - - + + • • • + + R n , 

\x x + 1 x + n-1/ 

mais 

logn = log - + log - + • • • + log — -, 

par suite 

12 1 3 1 4 

O(x) = — + log + log + log , 

x & 1 x-1 & 2 x + 2 & 3 

12 q- 1 q 

et, en y remplaçant x successivement par les fractions -, -,..., , - 

q q q q 

on tire de cette égalité 

^/1\ 1 , 2 q 3 q â 

%) =~î +l °H-q-ri +l ^2-2q-+-l + l0 H--' 

^/2\ <7 .2^.3 q 4 

%) = ~2 +l0 H-q-+-2 +l °Z2-2q-+-2 + l0 H--' 



O = h log h log h log - 

\ q j q-1 h l 2q-l h 2 3q - 1 5 3 

Jl) = _ï + l 0g 2 --± + log|-f + log^ 

\<7/ q b l 2q b 2 3q & 3 
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D'autre part, si l'on pose 

271 

co = — , 

en représentant par cp l'un quelconque des angles 

co, 2co, . . . , (cj - l)co, 

on vérifie sans peine les relations 

cos cp - cos(q + l)cp - cos(2q + l)cp - . . . , 
cos 2cp - cos(q + 2)cp - cos(2^ + 2)cp = ..., 



1 - cos qcp - cos 2qcp -cos3qcp -■■■, 

ainsi que la suivante : 

cos cp + cos 2cp + ■ ■■ + cos qcp - 0. 

On a donc 

/l\ q 2 

cos cp O - - - - cos cp + cos cp log - 

q 3 
cos(q + \)cp + cos cp log - - ■■■ 

(2\ q ,2 

cos 2cp O - - - - cos 2cp + cos 2cp log - 



,<?/ 



q . 3 

cos(^ + 2)<p + cos 2<p log - 



q + 2 Xl /r r °2 



/fl-l\ fl 2 

cos(fl - \)cp O = cos(fl - \)cp + cos(q - l)cp log - 



q t 3 

cos(2^ - l)cp + cos(q - l)cp log , 



2q-l v ' /r v/ /r °2 

? 



cos ^<p O - - - - cos qcp + cos ^<p log - 



- — cos 2qcp + cos qcp log , 



d'où, par addition 



cos cp O - J + cos 2cp O - ) + • • • + cos qcp O - 

- -q (cos cp + - cos 2cp + - cos 3cp + ■■■). 
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Il s'agit maintenant de chercher une expression différente de la somme 
de la série trigonométrique convergente qui constitue le second membre ; 
or, de l'identité 

1 - x cos <p - x n cos n<p + x n cos(n - l)cp 
- (1 - 2x cos (p + x ) 1 + x cos <p + x cos 2<p H + x M cos(n - l)<p , 

on déduit, pour toutes les valeurs de x ayant un module moindre que 
l'unité, 

l-xcos<» y o 

- - 1 + x cos cp + x cos 2<p + x cos 3cp H , 

1 - 2xcos<p + x z 

c'est-à-dire 

x-coscp ? 

- - cos <p - x cos 2cp - x cos 3cp — ••• . 

\-2x cos (p + x z 

Le premier membre est la dérivée de la fonction 

y - log a/1 - 2x cos <p + x 2 
qui s'annule avec x; d'ailleurs, pour x - 1, la série entière 

2 3 

-- cos (p — — cos 2cp — — cos 3cp 

est convergente, et si l'on égale sa somme calculée précédemment à la 
valeur de y pour x - 1, on obtient, d'après le théorème d'Abel, 

cos (p O - + cos 2(p O - + • • • + cos qcp O - - - log(2 - 2 cos cp). 

Si l'on considère à présent la somme 

cos pco cos rco + cos 2pco cos 2rco H h cos qpco cos qrco, 

où r est l'un quelconque des nombres 1, 2, 3, . . . ,q, on constate que cette 
somme est nulle pour toutes les valeurs de r autres que p et q - p, ces 

dernières valeurs la rendant égale à - ; en effet, cette somme a pour 
expression 

p + r . „ N p + ?" . p-r . „ N p-?" 

-^ smq——co cos(q + l)——co -y smq-—co cos(q + l)——co 

9 .p+r 9 . p—r 

*• sin — —co *• sm co 

2 2 
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Alors, si, dans la relation finale qui a été établie entre les nombres O I - 

O I -),... , O I - 1, on donne à (p successivement les valeurs co, 2co, ... , 

(q - l)co, en multipliant les équations formées de cette manière par cos pco, 
cos 2pco, ... , cos(q - l)pco, et les ajoutant à la relation 

»(i) + »(?) + ... + «(î) = ,| M(1) _, l0 g, / 

on trouve 

Ol- — -) + o(- = 20(l)-21og<7 + cospû; log(2 - 2 cos a») 

+ cos 2pco log(2 - 2 cos 2<y) H 

+ cos(q - l)pco log 2-2 cos(*7 - l)co ; 

les deux premiers termes étant mis à part, les coefficients des logarithmes 
également distants des extrêmes sont deux à deux égaux et, dans le cas 
où q est pair, le terme du milieu est 



cos 



qco ( qco\ 

Py log (2 -2 cos —y 



il a pour valeur 2 log 2 ou -2 log 2, suivant que p est pair ou impair. 
Enfin, d'après la relation des compléments, 

<D(ÎZZ|_<D|H = 7zcot-7z; 
\ <7 / V\) c\ 

on en conclut que, pour une valeur impaire de q, on a 

(p\ n P P / 2n\ 

O - =0(1)- -cot-7i-logfl + cos2-7i log 2-2cos — +■■■ 
\q) 2 q hH q h \ q) 

p I q-1 \ 

+ cos(q - l)-7i log 2 - 2 cos n , 

1 \ il 

et pour une valeur paire 

O - =0(1)- - cot-7i-logfl-cos2-7i log 2-2cos— +•• 
\q) 2 q h ^ q h \ q) 

+ cos(q-2)-n log 2-2cos- n\ + (-lf log 2. 

1 \ 1 I 
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Telles sont les formules qui permettent d'exprimer O(x) sous forme 
finie pour toute valeur rationnelle de x non égale à un entier négatif. 

Étude de l'équation O(x) - 0. — Si l'on se reporte au développement 
de O(x), on voit que le coefficient de i dans 0(a + Bi) a pour expression 



Il — *-*V 



-i (a + n) 2 + B 2 

n=0 v ' r 

Il ne peut être nul que pour 6 - ; la fonction O(x) n'admet donc pas de 
racines imaginaires. 
On a, d'autre part 

O(x) - O(x ) 



= 1 



x - xq *—* (x + u)(xq + n) 

Les termes du second membre sont tous positifs si les nombres x et Xq 
appartiennent simultanément à l'un quelconque des intervalles 

(0,+co), (-1,0), (-2,-1), ..., 

et comme 

0(+oo) - +oo, 0(-n - 0) - +oo, 0(-n + 0) - -oo, 

on eu conclut que, dans chacun des intervalles considérés, quand la 
variable x croît de la limite inférieure à la limite supérieure, la fonction 
O(x) augmente constamment de -oo à +oo et, par suite, s'annule une fois 
et une seule. Ainsi O(x) possède une racine positive et une infinité de 

racines négatives comprises respectivement entre et -1, -1 et -2, 

La racine positive appartient à l'intervalle (1, 2), car 

0(l) = -p, 0(2) = l-p, 

résultats dont le premier est négatif, et le second positif. Legendre i 1 ) et 
Gauss ( 2 ) ont calculé la valeur de cette racine qui, limitée à sept décimales, 
est 

x- 1.4616321..., 

alors 

T(x) = 0.8856024...; 



(^Mémoires de la classe des Sciences mathématiques et physiques de l'Institut de France, 
année 1809, p. 490. — Traité des fonctions elliptiques, t. II, p. 435-436. 
( ) Cari Friedrich Gauss Werke, t. III, p. 147. 
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c'est le seul minimum de la fonction T(x) dans l'intervalle (0, +00). 

Soit maintenant —n + h la racine négative de O(x) comprise dans 
l'intervalle (-n, -n + 1). Si, dans la relation 



O(x) - 0(1 - X) - -71 COt TIX, 




on fait x - -n + h, on obtient 


0(n + 1 - h) - 7i cot 7i/i; 


mais 

11 1 

0(n + 1 - /?) - 0(1 - h) = - — - + - — - + • • • + -, 

\-h 2-h n-h 


et comme la différence 


0(1 -h)- 


r / 1 1 IV 

logn - - — - + - — - + • • • + 

. & \l-h 2-h n-h). 





s'annule pour n - 00, on peut poser 

0(n + 1 - h) - logn - e, 

le nombre e tendant vers zéro quand n devient infini ; ainsi 

, 1 n 

h — — arc tang 



n e + logn' 

ou bien approximativement, si n est assez grand, 

logn 
dans cette hypothèse, la racine considérée ayant pour expression 

1 



x - -n + 



logn' 



il en résulte que les zéros se rapprochent de plus en plus des pôles à 
mesure que ceux-ci s'éloignent de l'origine. C'est à Hermite ( x ) que l'on 
doit ce résultat. Bourguet ( 2 ) a également étudié la même question, et il a 
calculé les racines correspondant aux deux premiers intervalles ; en outre, 
il a montré que, si Xp et x„ + i sont deux racines relatives à deux intervalles 
consécutifs, la racine relative à l'intervalle suivant est comprise entre x„ + i 
et 2x p+ i - x p . 



() Journal fur die reine und angewandte Mathematik, t. XC, 1881, p. 336-338. 
( 2 ) Atti délia R. Accademia délie Scienze di Torino, t. XVI, 1881, p. 738-772. 
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Application. — Courbe figurative de la fonction gamma. — Soit 

y = T(x), 

la variable étant supposée réelle. 

Lorsque x varie de à +00 la fonction T(x) varie de +00 à +00 sans 
s'annuler et, pour une valeur de la variable comprise entre 1 et 2, passe 
par un minimum qui vient d'être déterminé. 

Lorsque x varie de à —00, si l'on considère les intervalles successifs 
(-1, 0), (-2, -1), (-3, -2), . . . , la fonction T(x) est alternativement néga- 
tive et positive dans chacun de ces intervalles et devient infinie aux 
extrémités. 

Au moyen des tables, et en se servant de la relation 



T(-x) = — T(l-x) 

x 

pour les valeurs négatives de x, il est facile de construire la courbe par 
points ; on obtient ainsi la figure suivante : 
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A mesure que x s'écarte de l'origine dans le sens négatif, les points 
de la courbe correspondant à des maximums et à des minimums se rap- 
prochent de plus en plus de l'axe horizontal et de l'asymptote verticale 
la plus éloignée de l'origine. 
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VI. 



DEVELOPPEMENTS EN SERIES ENTIERES. 



Développement de log T(l + x) en série entière. — On a vu que 



logT(l + x) - -px + 



£- log(l + M _ log ( 1 + £) 



+ • 



Pour toute valeur de x de module inférieur à l'unité, les termes de cette 
série, à partir du second, sont développables en séries entières absolu- 
ment convergentes 



2 1 °> 1 4 

x - - x + - x - 
3 4 



ï--H)-| 

x i /-, x \ l/ x \ 2 l/ x \ 3 l/*\ 4 
2-log(l + 2)=2(2) -3(5) +5(2) — -, 

/ 

X . / X\ 1/X\2 1/X\3 1/X\4 

^- io s i+ ^ =9U) -*u + zU -•••' 



d'ailleurs, si r est le module de x, la série positive de terme général 

1 /r \2 



1 / r \ 3 1 / r \ 4 



1 / r y \ tr 
2\n) 3\n 
est convergente et a pour somme 

logT(l-r)-pr; 

par suite, si l'on pose 

1 1 1 

Sm- 1 1 h 

" ^n 2 M 3 n 
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en ordonnant par rapport aux puissances croissantes de x, on trouve 

2 3 

logr(l + X) = -p^ + S 2 y - S 3 j + ■ ■ ■ , |Z| < 1. 

Ce développement était connu d'Euler ( x ) ; il reste convergent pour 
x = 1, les valeurs absolues de ses termes décroissant alors constamment 
et tendant vers zéro, par conséquent, d'après le théorème d'Abel, 

s 2 s 3 s 4 

^y-y + T-'"' 

formule à laquelle Euler était parvenu d'une tout autre manière ( 2 ). 

Les sommes S 2 , S3, . . . , qui figurent dans les coefficients du dévelop- 
pement de log T(l + x), ont été calculées avec seize décimales, d'abord 
par Euler ( 3 ) jusqu'à S^, puis par Legendre ( 4 ) jusqu'à S35, et enfin avec 
trente-deux décimales par Stieltjes ( 5 ) jusqu'à S7Q. 

Le développement de log T(l + x) en série entière n'est pas d'un usage 
commode pour les calculs numériques ; Legendre ( 6 ), afin de le rendre 
plus convergent, l'a transformé de la manière suivante : 

On a 

2 3 

logr(l + X) = -p± + S 2 y - S 3 j + ■ ■ ■ , 

2 3 

X X z X° 



log(l + x) = -- — + 



d'où 



2 3 

x ,„ „ .x° 



log r(l + x) = - log(l + x) + (1 - p)- - (1 - S 2 )- + (1 - S 3 )- 

de même 

2 3 

I0 g r(l-X)--I0g(l-X) + (1-P)^-(1-S 2 )y-(1-S 3 )y----, 



( 1 ) Institutiones Calculi differentialis, pars posterior, § 384, p. 800. 

( 2 ) Commentant Academix Scientiarum imperialis petropolitanse, t. VII, 1734-1735, p. 156- 
157. 

( ) Institutiones Calculi differentialis, p. 456. 

( ) Exercices de Calcul intégral sur divers ordres de transcendantes et sur les quadratures, 
t. II, p. 65. 

( 5 ) Acta mathematica, t. X, 1887, p. 300-302. 
( ) Traité des fonctions elliptiques, t. II, p. 433. 
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par suite 

io g r(i + X ) -io g r(i -x) = io g T -^ + 2(1 - P y- + 2(1 -s 3 y- 



1 - X _ ,_ X _ ,„ „ , X 3 



+ 



mais 



r(l + x) - xT(x), 

Ti 



r(x)r(i - x) 



smnx 



donc 



77X 

r(i + x)r(i - x) , 

Sin7TX 



de sorte que 



7TX 

logHl + x) + logHl - x) - log ; 

sin nx 

il en résulte 

1 77 Y 1 1 — Y T Y 

logr(l + x) = - log -^— + - log - + (1 - p)- + (1 - S 3 )^- + • • • . 

& 2 5 sinra 2 & 1+ï r 1 ô ' 3 

Cette série, dont Legendre a calculé les huit premiers coefficients avec 
douze décimales après avoir converti les logarithmes népériens en loga- 
rithmes vulgaires, est très convergente pour les petites valeurs de x ; en 

l 
particulier, si l'on y fait x - --, elle devient 

, . 3 1 S 3 -l 1 S5-I 
P- 1 - l °^2-^^~^~5 ' 

formule qui peut servir au calcul de p. Euler, en l'employant, a obtenu la 
valeur de cette constante avec douze décimales. Stieltjes ( l ), reprenant ce 
calcul, a obtenu p avec trente- trois décimales. 

On trouve dans les Exercices de Calcul intégral de Legendre les tables 
de log T(x) depuis 1 jusqu'à 2, de millième en millième, d'abord avec sept 
décimales (t. I, p. 302-306), puis avec douze décimales (t. II, p. 85-95). 
Gauss ( 2 ) a fait calculer par Nicolai les tables de log T(x) à vingt décimales, 
de centième en centième, pour les valeurs de x comprises entre 1 et 2. 



( x ) Acta mathematica, t. X, 1887, p. 302. 

( ) Cari Friedrich Gauss Werke. t. III, p. 161-162. 
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Enfin, Bellavitis ( x ) a donné les tables de logT(x) à huit décimales, de 
dixième en dixième, depuis 10 jusqu'à 12. 

1 
Développement de T(l + x) et de — en séries entières. — On a 

vu (p. 11) que la possibilité du développement de la fonction T(l + x) en 
série entière résulte de ce que son inverse est une fonction transcendante 
entière. Il ne nous reste donc plus qu'à montrer comment s'effectue le 
calcul des coefficients de chacune des deux séries en partant du déve- 
loppement de logT(l + x). 
Soit d'abord 

logT(l + x) = 1 + A x x + A 2 x 2 + ■■■; 

la fonction T(l + x) étant dérivable à l'intérieur de son cercle de conver- 
gence, on obtient en dérivant le développement de son logarithme 

Hl+x) „ „ 2 

c'est-à-dire 

A 1 + 2A 2 x + 3A 3 x 2 + ■ ■ ■ = (1 + A x x + A 2 x 2 + ■■■ )(-p + S 2 x - S 3 x 2 + ■■■), 

d'où 

M - -p, 

2A 2 = -pA x + S 2 , 

3A 3 = -pA 2 +A 1 S 2 -S 3 , 



nA n = -pA n _ x + A n _ 2 S 2 - A n _ 3 S 3 + ■■■ + {-ï) n S n , 



ces relations permettent de déterminer successivement les coefficients 

A\, A 2 , ..., A n , Les vingt premiers d'entre eux ont été calculés avec 

douze décimales par Jeffery ( 2 ). 
Soit maintenant 

1 +Bix + B 2 x 2 + •••; 



T(l + x) 



( ) Memorie del reale Istituto veneto di Scienze, Lettere ed Arti, vol. XVIII, 1874, p. 126-162. 
( ) The quarterly Journal ofpure and applied Mathematics, vol. VI, 1864, p. 96-97. 
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si l'on dérive les deux membres de cette relation, elle donne 



( r A ' =B 1 + 2B 2 x + 3B 3 x 2 + 



T(l+x) T(l+x) 

c'est-à-dire 

(1 + B x x + B 2 x 2 + ■ ■ ■ )(p - S 2 x + S 3 x 2 ) = B 1 + 2B 2 x + 3B 3 x 2 + ■■■ , 

et, en identifiant les coefficients des mêmes puissances de x dans les deux 
membres, on trouve 

2B 2 = P B 1 - S 2 , 

3B 3 = pB 2 -B 1 S 2 + S 3 , 



nB n - pB n _i - B n _ 2 S 2 + B n _ 3 S 3 (-1) S n , 



Ces équations se déduisent de celles qui ont servi à déterminer Aj, A 2/ 
A 3 , ... en y changeant les signes de p, S 2 , S 3/ ... ; elles donnent suc- 
cessivement les coefficients B\, B 2 , B 3 , Les vingt-deux premiers de 

ces nombres ont été calculés avec seize décimales par Bourguet ( x ). Il y 
a intérêt, pour simplifier les opérations, à mettre en évidence dans le 
développement le facteur 1 + x, ce qui revient à retrancher l'unité des 

coefficients p, S 2/ S 3/ On forme ensuite les neuf premiers multiples 

de ces nombres, puis, afin de vérifier les multiplications par le renver- 
sement des facteurs, on cherche également les neuf premiers multiples 
de chacun des coefficients B\, B 2 , B 3 , ... à mesure qu'on les obtient. On 
emploie la méthode abrégée pour effectuer les multiplications ( 2 ). 
On tire de l'identité 

1 = (1 + A x x + A 2 x 2 + ■ ■ • )(1 + B x x + B 2 x 2 + ■■■) 



les relations suivantes 



A 1 + B x = 0, 
A 2 +A 1 B 1 +B 2 = 0, 

A n +A n _iBi + --- + B n = Q, 



( 1 ) Annales scientifiques de l'École normale supérieure, 2 e série, t. X, 1881, p. 224. 
( ) Pour le détail des calculs on se reportera à la Thèse de Bourguet {Annales scienti- 
fiques de l'École normale supérieure, 2 e série, t. X, 1881, p. 209-212). 
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on peut ainsi obtenir les valeurs des coefficients B\, £>2, ... au moyen de 
celles des coefficients A\, Ai, . . . , ou inversement. 

La méthode que nous venons d'exposer pour établir les équations 
relatives aux coefficients des développements de la fonction T(l + x) 
et de son inverse diffère complètement de celle indiquée par Bourguet 
dans sa Thèse; elle constitue, croyons-nous, un réel progrès sur celle 
dernière ( a ). 

Développement de P(l + x) et de Q(l + x) en séries entières. — Les 

termes de la série 

P(l + x) = + 



1+x 1 2 + x 1-2 3 + x 

1 •2---(n - 1) n + x 
pour toute valeur de x vérifiant l'inégalité 

|*| < 1, 

sont développables en séries entières absolument convergentes 

= 1 -X +X -x + • • • , 



1 + x 



1 1 X X 2 X 3 






2 + x 2 2 2 2 3 2 4 



1 X X 2 X 3 



n + x n n 2 n 3 n 4 



d'ailleurs, si r est le module de x, la série positive de terme général 
1 1 1 /l r r 2 \ 

+ ^r + -=■ + ■■■ 



l-2---(n-l)n-r 1 • 2---(n - 1) \n n 2 n 3 
est convergente ; par suite, en posant 



(~lf 



1 1 1 

1- r + 



1 2P +1 1 • 2 3P+1 

y ' l-2---(n-l) n p+l 



( ) Voir une Note que nous avons publiée dans les Annales de la Vacuité des Sciences 
de Marseille, t. XI, 1901, p. 117-124. — Voir aussi la Thèse de Lindhagen, Studier ôjver 
gamma-funktionen, p. 12-13 et p. 19. 
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quand on ordonne par rapport aux puissances croissantes de x, on trouve 

P(l + x) - dQ + a\X + a 2 x H , 

série dont le rayon de convergence est l'unité (p. 21). Pour calculer les 

il 
coefficients a§, a\, ai, • • • , on commence par former les quotients -, — — , 

1 

, ... ; ensuite on les divise respectivement par 2, 3, 4, ... ; puis les 

résultats obtenus sont divisés de nouveau par 2, 3, 4, . . . , et ainsi de suite. 
Il est maintenant facile de déterminer le développement de la fonction 
Q(l + x). En effet, 

r(l + x) = 1 + A x x + A 2 x 2 + ■■■ , 
P(l + x) - a$ + a^x + a 2 x + ■ ■ ■ , 



par conséquent, 

Q(l + x) - 1 - a + (Ai - a{)x + (A 2 - a 2 )x 2 + • • • ; 

tel est le développement de Q(l + x) en série entière. On en déduit sans 
peine le développement de la fonction Q(x) elle-même, car 

Q(l + x) - xQ(x) + -, 
e 

d'ailleurs 

«o = l — / 
e 



d'où 



Q(x) =A 1 -a 1 + (A 2 - a 2 )x + (A 3 - a 3 )x 2 + 



développement qui, comme on le sait (p. 21), converge pour toute valeur 
dex. 

Les coefficients A\-a.\,A 2 -a 2 ,..., jusqu'à A 18 - a 18 / on t été calculés 
par Bourguet avec seize décimales ( a ). 

Développement de 0(1 + x) et de ^(l + x) en séries entières. — 

Les fonctions 0(1 + x) et ^(l + x) sont développables en séries entières 
de rayon de convergence égal à l'unité; en effet, par deux dérivations 
successives du développement 

2 3 

io g r(i + x) - - P - + s 2 — - s 3 — + • • • , |*| < î, 



( ) Annales scientifiques de l'École normale supérieure, 2 e série, t. X, 1881, p. 227. 



— 69 — 

on obtient 

0(1 + x) = -p + S 2 x - S 3 x 2 + ..., \x\<l, 
W(l + x) = S 2 -2S 3 x + 3S 4 x 2 - ..., \x\<l. 

Gauss ( x ), dans son Mémoire sur la série hypergéométrique, a donné, 
d'après les calculs de Nicolai, en même temps que les valeurs de log T(x), 
celles de O(x) de 1 jusqu'à 2, de centième en centième, et avec dix-huit 
décimales. 



( 1 ) Cari Friedrich Gauss Werke, t. III, p. 161-162. 
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VIL 

APPLICATIONS. 



M. Appell a signalé de très intéressantes applications des fonctions 
T(x) et W(x) à l'évaluation de la limite de produits ou de séries dont le 
terme général est une fonction rationnelle de l'indice, et à la résolution de 
certaines équations fonctionnelles, où figure une fonction rationnelle de 
la variable. On est ainsi conduit à des résultats qui sont susceptibles d'im- 
portantes généralisations. C'est ce dernier problème que M. Appell ( x ) a 
étudié à l'aide de fonctions nouvelles. 

Nous nous proposons, dans le présent Chapitre, d'exposer en détail 
les applications indiquées par M. Appell, et nous terminerons en donnant 
la résolution de l'équation de Lindhagen et des équations de Crelle, 
questions qui appartiennent au même ordre d'idées. 

Limite de produits infinis convergents dont le terme général est 
une fonction rationnelle de l'indice. — Soit 

(n + a) a {n + bf ■■■(n + l) A 

1j — A 

{n + a') a '{n + b'f ■■■(n + l') A ' ' 

en désignant par A, a, b,..., I, a' ', b' ,... , V des constantes quelconques, et 
par a, j8, . . . , A, a! ', §',..., A' des exposants entiers et positifs. Si l'on pose 



P n - UqU\ ■ ■ ■ u 



h ii 



( ) Les recherches de M. Appell, publiées d'abord dans les Comptes rendus hebdoma- 
daires des séances de l'Académie des Sciences, t. LXXXVI, 1878, p. 953-956, et t. LXXXIX, 
1879, p. 841-844 et p. 1031-1032, ont été réunies par lui en un Mémoire inséré dans les 
Mathematische Annalen, t. XIX, 1882, p. 84-102. A ce propos il faut citer un travail de 
Scheibner, curieux et peu connu, qui parut en 1861 dans les Berichte iiber die Verhand- 
lungen der kôniglich sachsischen Gesellschaft der Wissenschaften zu Leipzig. Mathematisch- 
physische Classe, t. XIII, p. 64-135. 
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on a 



\a (a + 1) •• • (a + n - 1) 
P„ - A nl - i — 

[fl / (fl / + l)---(fl / + n-l)] a 
[b (b +l)---(6 +n-l)f •••[/(/ +l)---(/ +n-l)] A 

[&'(&' + 1) •■•(&'+ n - l)f ' • • • [/'(/' + 1) •••(/' + n - 1)] A ' 

or 

T(x + n) = x(x + 1) • • • (x + n - l)r(x), 

par suite 

T a '(a')rP'(b')---T A '(l') T a (a +n)TP (b +n)---T A (l +n) 



P n — A 



T a (a)TP (b)---T A (l) Y a \a' + n)lP'(b' + n)---T A '(l' +n)' 



telle est l'expression du produit des n premiers facteurs. Il faut, pour 
que le produit P n soit convergent, que son terme général u n tende vers 
l'unité. On doit donc avoir 



et 



Soient alors 



A = l, 



a + ]3 + --- + A = tt / + J S / + --- + A / . 



s - a a +b f> H vl À, 

s' - a' a' + b'p' + ■■■ + l'A'; 



on sait que, pour n - oo, le rapport 

1 T(x + n) 



n x i . 2 . . . ( M _ 1) 



devient égal à l'unité ; en tenant compte de ce résultat, on constate, après 
une transformation facile, que le produit P n ne peut être convergent que 
pour 



s - s', 



et, cette condition étant réalisée, la limite P de P n est donnée par la 
formule 

T a \a')TP\b')---T A \r) 



P = 



T a (a )T? (b )---T A (/)" 
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C'est à de Gasparis ( x ) qu'est due la première étude des produits dont 
le terme général est une fonction rationnelle de l'indice. Mais, comme l'a 
remarqué Enestrôm ( 2 ), les conditions pour qu'un produit de ce genre ait 
une limite étaient connues de Stirling( 3 ). 

On déduit, sans difficulté, des résultats précédents les conditions 
relatives à la convergence absolue d'une série dans laquelle le rapport 
d'un terme au précédent est une fraction rationnelle telle que 

u n+ i (n + a ) a (n + b )P ■ ■ ■ (n + l ) A 



u n (n + a') a ' (n + b'f ■ ■ ■ (n + l') A " 

le numérateur et le dénominateur ayant le même degré. En effet, si l'on 

pose 

_ U\ u 2 u n 
In— ' 



on a 



Soient alors 



Uq ll\ U n —\ 

u n - UqP h . 

o + t i - a a +b fi +••• + / À, 
a' + %'i = a' a' + b'p' + ■■■ + l'A'; 

d'après ce que l'on vient de voir, le module de u n peut se mettre sous la 
forme 

\Un\ — AnH , 

en désignant par A n un coefficient positif qui tend vers une limite non 
nulle pour n - oo. Ainsi 

y,0'-0 I,. I _ 1 

ri \Un\ — <\-nr 

et l'on en conclut que la série de terme général u n n'est absolument 
convergente que pour 

o' - o > 1. 

Si le rapport d'un terme au précédent a pour expression 

u - nP + anP~ l + ainP~ 2 + ■ ■ ■ + a p _i 
u n nP + bnP~ l + binP~ 2 + • • • + b v _\ 



( 1 ) Rendiconto délie adunanze e de' lavori délia reale Accademia délie Scienze di Napoli, 
2 e série, 1. 1, 1852, p. 199-201. 

( 2 ) Ofversigt afkongl. Vetenskaps-Akademiens Fôrhandlingar, t. XXXVI, 1879, n° 9, p. 71- 
84. 

(^)Methodus differentialis : sive Tractatus de summatione et interpolatione serierum infini- 
tarum, p. 37-40. 
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la condition de convergence absolue est que la partie réelle de la diffé- 
rence b - a soit supérieure à l'unité. Cette règle de convergence, qui est la 
généralisation de celle donnée par Gauss pour les séries à termes positifs, 
a été utilisée par Hermite ( x ) dans certaines applications du théorème de 
Mittag-Leffler. 

Sommation des séries dont le terme général est une fonction ration- 
nelle de l'indice. — Si l'on considère une série dont le terme général u n , 
fonction rationnelle de n, puisse être décomposé en fractions simples de 
telle sorte que 



Un 



-L 



Ai A 2 

+ 



n + a (n + a) 2 



+ ... + 



A, 



(n + a) c 



en posant 

Sfi — Uq + U-y + • • • + U n —\, 

on voit que S n est la somme des expressions 



L 



Ai 



1 1 

+ 



a a + 1 



+ ••• + 



la 2 (a + l) 2 



1 




a + n-1 
1 


/ 


(a + n-1) 2 



M* 



+ 



(a + iy 



-+••• + 



Or 



(a + n- l) c 



( 1)W \^ n \x + m)-^ n \x)] = -!— + — 

1 -2---n L J x n+1 (x + l) n 

par suite 



+i ' (x + tn-l) n+1 ' 



1 
Sn - EAx[o(fl + n) - 0(a)] - -£A 2 [o> + n) - ®'(a)] + 



+ 



(-1) 



a—\ 



YA a \^- l \a + n) - O^" 1 )^)]; 



l-2---(a-l) 

telle est l'expression de la somme des n premiers termes. Il reste à cher- 
cher sa limite, pour n - oo, lorsque la série est convergente. 



(- 1 ) Faculté des Sciences de Paris. — Cours de M. Hermite rédigé en 1882 par M. Andoyer, 
3 e éd., p. 136-137. 
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Si l'on réduit les fractions simples au même dénominateur l'expres- 
sion de u n devient 

a nP + ainP~ l + ■■■ + aP 
Un ~ \ • 

b n1 + b\rfi i + • • • + W 

On doit écarter, d'abord, l'hypothèse p - q ^ 0, car, dans ce cas, u n ne 
tend pas vers zéro quand n augmente indéfiniment ; si l'on suppose alors 
p - q < 0, on peut mettre u n sous la forme 

1 

^ a + -(ai + a n ) 

b 0+-(h+pn) 

les nombres a n et § n s'annulant pour n - oo ; il en résulte que, à partir 
d'une valeur suffisamment grande de n, on aura 



a 

u n 



b nl-P 



< 



n q-p+l' 



en désignant par k une constante positive. La série de terme général u n 
est donc encore divergente pour q - p - 1 ; elle ne converge que si l'on 
a p - q ^ -2. Cette dernière condition étant remplie, les coefficients tels 
que Ai vérifient la relation 

LAi = 0, 

puisque, lorsqu'on réduit les fractions simples au même dénominateur, 
le premier membre de l'égalité précédente est le coefficient de rfl dans 
le polynôme du numérateur; par conséquent, si n croît indéfiniment, 
l'expression 

XA^OCa + n) - O(fl)] - LA 2 [o(a + n) - O(a) - logn] 

a pour limite 

-ZAi®{a); 

on en déduit que la série considérée a pour somme 

i (—l) a 

S = -EAiO(a) + -ZA 2 0'(a) - • • • + K / -ZA a ^ a - l \a). 

1 1 • 2- •• (a - 1) 

Dans le cas particulier où le dénominateur de la fraction rationnelle 
u n n'admet que des racines simples, cette formule se réduit à 

S = -EAiO(a), 
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ou bien 



S = -£Ai[<D(a)-0(l)]. 



Ainsi, lorsque toutes les racines sont rationnelles, la somme de la sé- 
rie est réductible aux logarithmes et aux fonctions circulaires. C'est à 
M. Appell i 1 ) que l'on doit cette remarque. 

Eduard Weyr( 2 ), dans un Mémoire sur la fonction O(x), a donné 
plusieurs exemples numériques à propos de la méthode de sommation 
que nous venons d'exposer. Nous citerons seulement, d'après Jensen ( 3 ), 
le développement suivant 

n=co „ r 

O(a) 



y 1 

*—* (n + a)(n + &)■■■ (n + /) 



(b - a){c -«)■••(/-«) 



0(b) 0(1) 



(a - b)(c -b)---(l-b) (a- l)(b -/)••• (Jfc - /) 

La sommation d'une série entière, où le coefficient de x n est une 
fonction rationnelle de n, est un problème plus général, mais moins 
simple, que celui qui vient d'être traité. Il fait l'objet de l'un des Mémoires 
posthumes d'Abel ( 4 ). Hermite ( 5 ) s'est aussi occupé du même sujet, dans 
le cas où la fonction rationnelle se réduit à un polynôme. 

Résolution de certaines équations fonctionnelles. — Nous nous 
bornerons à résoudre deux questions qui présentent la plus grande ana- 
logie avec celles précédemment étudiées. 

I. Soit d'abord l'équation fonctionnelle 

F(x + 1) = R(x)F(x), 

où R(x) désigne une fonction rationnelle de x. On peut mettre cette fonc- 
tion R(x) sous la forme 

(x + ar(x + bf---(x + l) A 
K(X) — A- 



(x + a') a ' (x + b'f ■ ■ ■ (x + V) A " 



en désignant par A,a,b,.. . ,l,a' ,b' ,. . . ,V des nombres indépendants de x, 
et par a, §,..., A, a' , §',..., A' des exposants entiers et positifs. 



( ) Comptes rendus hebdomadaires des séances de l'Académie des Sciences, t. LXXXVI, 1878, 
p. 954. 

( ) Casopis pro péstovâni mathematiky afysiky, t. XXI, 1892, p. 176-180. 

( 3 )Nyt Tidsskrift for Mathematik, t. II B, 1891, p. 49. 

( ) Œuvres complètes, éd. L. Sylow et S. Lie, t. II, p. 14-18. 

( ) Bulletin de la Société physico-mathématique de Kasan, 2 e série, t. VIII, 1899, p. 107-108. 
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Or, on a successivement 

F(x + 1) = £(x)F(x), 

F(x + 2) = R(x + l)F(x + 1), 

F(x + n) = R(x + n- l)F(x + n-ï), 

par suite 

F(x + n) = R(x)R(x + l)---R(x + n- l)F(x), 

et, comme on l'a vu déjà, on trouve sans peine 

„ T a ' (x + a')!?' (x + b')---T A '(x + V) 



F(x + n)=A r 



T a (x + a)TP(x + b)---T A (x + l) 



T a (x + a + n)TP(x + b + n)---T A (x + l + n) 
T a '(x + a' + n)lF{x + V + n) ■ ■ ■ T A '(x + V + n) 

Une solution particulière de cette équation est évidemment 

= A% r*(x + a)r% + fr)---r A (x + o 

r«'(x + a')TF(x + ¥)■■■ T A '(x + V) 
et la solution générale sera 

F(x) - 0(x)S(x), 

en représentant par 0(x) une fonction périodique arbitraire de période 
égale à l'unité. 
Soient 

s = aa + bp + ■■■ + lA,s' = a' a' + b'fi + ■■■ + V 'A' ', 
p - a + 18 + • • • + A, p' = a' + fi' + ■ • ■ + A'; 

si l'on impose à la fonction F(x), non seulement la condition de vérifier 
l'équation fonctionnelle donnée, mais encore celle d'être telle que le 
rapport 

1 1 F(x + n) 

Â» ' n (p-p')x+SS' ' r 2 ~ _ Ip-P' 

ait une limite K pour n - oo, il n'y aura plus qu'une seule solution 

r a (x + a)r% + fr)---r A (x + /) 



F(x) = K 



T a ' (x + a')TF (x + &') • • • P 1 ' (x + V) ' 
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Mellin ( 1 ), dans un Mémoire sur les séries hypergéométriques d'ordre 
supérieur, a donné le nom de fonction gamma à toute solution de l'équation 
fonctionnelle 

F(x + 1) - A (* + <0(* + ft)-"(* + Q r(x) 

c'est-à-dire à toute fonction de la forme 

_ T(x + a)T(x + b)...T(x + l) 
W T(x + a')T{x + ¥)■■■ T(x + /') 

la fonction &(x) admettant l'unité comme période. En se plaçant à ce 
point de vue très général, toute fonction rationnelle est une fonction 
gamma. 

Une équation fonctionnelle semblable à celle qui vient d'être étudiée 
est la suivante 

F(x + 1) = r(x)F(x), 

où r(x) est une fonction de la forme 

, x m (x + af(x + bf---{x + l) A 
r(x) - e 



{x + a') a \x + b'f ■■■{x + l') A " 

si &(x) est une fonction ayant l'unité pour période, la solution générale 
de cette équation est 

F(x) - 0(x)S(x), 

en posant 

T a (x + a)TP(x + b)---T A (x + l) 



S(x) - e 



MX 



T a '(x + a')TP\x + &')••• tA '( x + l ')' 



fonction qui, comme nous l'avons déjà dit (p. 21), a été de la part de 
Mellin ( 2 ) l'objet d'une étude très étendue. 

II. Soit en second lieu l'équation fonctionnelle 

F(x + 1)= F(x) + R(x), 

dans laquelle R(x) désigne une fonction rationnelle donnée susceptible 
d'être mise sous la forme 



R(x) - Yu 



A 1 A 2 A c 

+ ^7 + • • • + 



x + a (x + a) 2 (x + a) c 



( ) Acta Societatis Scientiarum fennicx, t. XXIII, 1897, n° 7. 
( 2 ) Acta mathematica, t. VIII, p. 37-80. 
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On a successivement 

F(x + 1) = F(x) + R(x), 

F(x + 2) = F(x + 1) + R(x + 1), 

F(x + n) - F(x + n - 1) + R(x + n - 1), 

par suite 

F(x + n) - F(x) - K(x) + R(x + 1) + • • • + R(x + n - 1), 

et, en raisonnant comme on l'a fait dans le cas d'une série, on trouve 

F(x + n) - F(x) - L^ijOC* + a + n) - 0(x + a)l 

Z^lo^x + a + n) - 0'(x + a)]+ • • • 

(_-\ \a-l 
+ 1 . 1 .\ a _ 1) ^ A <¥ {a ~ l) ( x + «"«)- <^ a ' 1 \x + a)}. 

On voit tout de suite qu'une solution particulière est 
S(x) - 2>!0(x + a) 

1 / 1 \ct— 1 

- tLA 2 &(x + «) + ■■■ + — ^ rrLAaO^-^Cx + a), 

1 1 • 2- ••(a - 1) 

de sorte que la solution générale est donnée par la formule 

F(x) - 0(x) + S(x), 

en désignant encore par 0(x) une fonction périodique de période égale 
à l'unité. 

Si l'on suppose le degré du dénominateur de la fraction rationnelle 
R(x) supérieur de deux unités au moins au degré du numérateur, et 
que l'on assujettisse la fonction F(x), non seulement à vérifier l'équation 
fonctionnelle, mais, en outre à être telle que F(x + n) tende vers une limite 
K pour n - oo, il n'y a plus, dans ce cas, que la solution unique 

F(x) = K + 2>!0(x + a) 

- tLA 2 ®'(x + «) + •■■ + — K —t- -LA a ^ {a - l \x + a). 

1 1 • 2- •• (a - 1) 
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Équation de Lindhagen. — L'équation fonctionnelle 

F(x + 1) = xF(x) + R(x), 

où R(x) désigne un polynôme entier en x, se rattache par sa forme aux 
équations dont nous venons de nous occuper. Cette équation a été étudiée 
par Lindhagen ( x ), mais c'est Jensen ( 2 ) qui en a fait connaître la solution 
générale sous forme explicite, et voici de quelle manière : 
De l'égalité 

F(x + 1) F(x) R(x) 



on tire 



d'où 



T(x + 1) T(x) T(x + 1)' 

F(x + 2) F(x + 1) _ R(x + 1) 
T(x + 2) ~ T(x + 1) " T(x + 2)' 

F(x + n) F(x + n - 1) R(x + n - 1) 
T(x + n) T(x + n - 1) T(x + n) 

F(x + n) F(x) v 1 " R(x + p - 1) 



E 



T(x + n) T(x) ^ T(x + p) 

Or, pour toute valeur de x non égale à un entier négatif, la série de 
terme général 

R(x + n - 1) 

T(x + n) 
est convergente, le rapport de ce terme au précédent s'annulant quand 
n devient infini ; d'ailleurs, si m est le degré du polynôme R(x), on peut 
poser 

R(x) - Aq+ A\X + A 2 x(x - 1) + • • • + A m x(x - 1) • • • (x - m + 1), 

car il suffit de donner à x les m + 1 valeurs 0, 1, 2, ... ,tn pour déterminer 
les coefficients Aq, Aj, A^, ■ ■ ■ , A m . Ainsi 

R(x) A | A 1 | A 2 | | A m 



r(x + i) r(x + i) r(x) r(x-i) r(x-m + i)' 

R(x + 1) A | Ax | A 2 | | A m 



r(x + 2) r(x + 2) r(x + i) r(x) r(x-m + 2)' 



( ) Studier ôfver gamma- funktionen och nâgra besïagtàde transcendenter, p. 44-55. 
( 2 )Nyt TidsskriftforMathematik, t. II B, 1891, p. 60-62. 
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par suite (p. 22) 



P(x) P(x - 1) P(x - m) _ 1 y. R(x + n - 1) 

°f(x) + 1 T(x-l) + '" + m T(x-m) " ë 2-j f(x + n) ' 



mais (p. 20) 



d'où 



eP(x) = -1 - (x - 1) - (x - l)(x - 2) 

- (x - Z)(x - 2) • • • (x - p + 1) 
+ (x - l)(x - 2) • • • (x - p) e P(x - p), 
T(x) - (x - p)(x - p + 1) • • • (x - 1) T(x - p), 



T(x)P(x-p) 
e—- -!- = eP(x) + X p , 

T(x -p) v 



en représentant par Xp le polynôme 

1 + (x - 1) + (x - l)(x - 2) + • • • + (x - l)(x - 2) • • • (x - p + 1); 
si donc on désigne par S(x) la somme de la série 

R(x + n - 1) 



I 



T(x + n) 



on a 



S(x)r(x) =A X X X + A 2 X 2 + A 3 X 3 + ■■■+ A m X m + e(A +A 1 + ---+A m ) P(x). 

Il est facile maintenant de résoudre l'équation fonctionnelle donnée. 
En effet, le rapport 

F(x + n) 

T(x + n) 

a nécessairement une limite pour n = oo alors, si l'on considère une fonc- 
tion G(x) qui vérifie l'équation fonctionnelle et soit telle que le rapport 

1 G(x + n) 



n x 1 • 2 • • • (n - 1) 



ait pour limite une constante K quand n croît indéfiniment, cette fonction 
G(x) aura pour expression 

G(x) - KT(x) - S(x)r(x); 
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d'autre part, des égalités 

F(x + 1) = xF(x) + R(x), 
G(x + 1) = xG(x) + R(x), 

on déduit 

F(x + 1) - G(x + 1) = x[F(x) - G(x)], 

et, si &(x) est une fonction admettant une période égale à l'unité, la 
solution générale a pour expression 

F(x) = G(x) + &(x)T(x), 

c'est-à-dire 

F(x) = [K + &(x)]T(x) - S(x)T(x), 

ou bien, en négligeant la constante K qui ne modifie en rien la forme de 
la fonction arbitraire &(x) 

F(x) = &(x)T(x)-(A 1 X 1 +A 2 X 2 + ---+A m X m )-e(A +A 1 + ---+A m )P(x). 



Équations de Crelle. — Soient x et y deux variables, et z un entier 
positif : si l'on pose 



F(x, y, z) - x(x + y){x + 2y) ■ ■ ■ (x + z - \y), 

en désignant par a un entier positif et par À un nombre quelconque, on 
voit immédiatement que la fonction F(x, y, z) vérifie les équations 

F(x, y,z + a) = F(x, y, z) F(x + zy, y, a), 

F(Ax, Ay, z) - A z F(x, y, z), 

F(x,y,l) = x. 

Lorsque l'on fait y - dans F{x,y,z), cette fonction se réduit à l'ex- 
ponentielle x z ; les relations précédentes correspondent donc aux trois 
propriétés caractéristiques de l'exponentielle x z exprimées par les égali- 
tés : 



x z+a 


- x z x\ 


(Ax) z 


= A z x z , 


x 1 


- x. 
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La fonction x z ayant une signification précise pour toute valeur de l'ex- 
posant z, on est amené par analogie à se demander s'il n'existe point une 
fonction F(x, y, z) qui, pour toute valeur entière et positive de z, soit re- 
présentée par le produit déjà considéré, et satisfasse aux trois équations 
qu'il vérifie. On a désigné, d'après Kramp, une telle fonction sous le nom 
de faculté analytique. L'étude des fonctions de ce genre, ébauchée par Van- 
dermonde if), a été l'objet de nombreux travaux de la part des analystes 
allemands de la première moitié du xix e siècle ; il faut citer, entre autres, 
les recherches de Kramp ( 2 ), Bessel( 3 ), Crelle, Oettinger, Ohm( 4 ). Mais 
les résultats obtenus étaient loin d'être satisfaisants lorsque Weierstrass, 
dans un Mémoire ( 5 ) dont nous avons eu plusieurs fois l'occasion de par- 
ler, parvint à élucider toutes les difficultés par une méthode aussi sûre 
qu'élégante. Voici quel a été son point de départ. 

Le problème dont il s'agit revient à chercher la solution générale des 
trois équations suivantes, dites équations de Crelle ( 6 ), 

F(x, y,z + a) - F(x, y, z) F(x + zy, y, a), 
F(ax, ay, z) - a z F(x, y, z), 
F(x, y,l) = x. 

Or, si l'on permute z et a dans la première équation, elle devient 

F(x, y,z + a) - F(x, y, a) F(x + ay, y, z), 

d'où, en remplaçant x par x - ay, 

F(x-ay,y,z + a) 

F(x,y,z) = — — . 

F(x - ay, y, a) 

Soit alors 

x- ay - ty; 

on trouve 

c , v F(ty,y,z + a) 

F(x,y,z) - 



F(ty,y,a) 



(^Mémoires de Mathématique et de Physique, tirés des registres de l'Académie royale des 
Sciences, 1772, l re Partie, p. 489-498. 

( 2 ) Analyse des réfractions astronomiques et terrestres, Strassbourg, 1799, in-4°. 

( 3 ) Abhandlungen von Friedrich Wilhelm Bessel, éd. Engelmann, t. II, p. 342. 

( ) Les Mémoires de Crelle, Oettinger et Ohm ont été publiés dans le Journal fur die 
reine und angewandte Mathematik, t. VII, 1831 ; t. XXXIII, 1846 ; t. XXXV, 1847 ; t. XXXVIII, 
1849; t. XXXIX, 1850; t. XLIV, 1852. 

( ) Ùber die Théorie der analytischen Facultaten (Mathematische Werke, 1. 1, p. 153-221). 

() Journal fiir die reine und angewandte Mathematik, t. VII, 1831, p. 269. 
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d'où, en tenant compte de la seconde équation, 



par suite, si l'on pose 



F(t,l,z + a) 



%)= F K4 



on en déduit 



F(x,y,z) = y 



.±3 



Inversement, il est facile de constater que, quelle que soit la forme 
de la fonction arbitraire G(x), toute fonction F(x, y, z) définie par cette 
relation satisfait aux deux premières équations de Crelle. En effet, 

G\-+z + a) G\-+z) G\^^- + a) 

.z+a \y i _* \v in \ v I 



F(x,y,z + a) = y = y --— y 



c'est-à-dire 
d'autre part 



ri I ' r( x \ r ( x + z y 

\y) \y) \ y i 



F(x, y,z + a) = F(x, y, z) F(x + zy, y, a); 



F(ax,ay,z) - (ay) 



AH 



^ 



c'est-à-dire 



F(ax, ay, z) = a z F(x, y, z). 



Maintenant, pour que la troisième équation de Crelle soit vérifiée, il 
faut que l'on ait 

y — t^ = x ' 
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ou bien 

lx \ _ x r (x\ 

\y + ! y \y)' 

équation dont la solution générale est 

G ÊK) r Ê)' 

en désignant par &(x) une fonction arbitraire de x admettant une période 
égale à l'unité ; par conséquent l'expression la plus générale des fonctions 
F(x, y, z) vérifiant les trois équations de Crelle est donnée par la formule 

r(* + z)e(* + z) 

r, \ z \y / \y / ,n 





( ) Weierstrass, Mathematische Werke, 1. 1, p. 162. 



TABLE METHODIQUE. 



i. 

HISTORIQUE. 



Pages. 

Wallis et Stirling 1 

Euler 1 

Gauss 2 

Legendre 3 

Période moderne 3 

Weierstrass 4 

Prym, Hermite 5 



IL 



ETUDE GENERALE DE LA FONCTION GAMMA. 

Définition de la fonction gamma 6 

Constante d'Euler 7 

Fonction factorielle de Weierstrass 8 

Formule de Weierstrass 8 

Fonctions de Heine, Barnes, Appell, Pincherle, Alexejewsky 9 

Holomorphie de l'inverse de T(x) 10 

Relation fonctionnelle 11 

Module de T(a + fii) 12 

Résidus de T(x) 13 

1 T(x + ri) 

Limite pour n - oo de l'expression — - — - — — 14 

Étude d'une série de Stirling 15 

Convergence de la série hypergéométrique 16 

Fonctions de Prym 17 



86 



Fonction de Mellin 21 

Problème de M. Appell 21 

Fonction de Bourguet 22 

Étude de l'équation P(x) = 23 



III. 



PROPRIETES DE LA FONCTION GAMMA. 

Relation des compléments 26 

Formule de Legendre 27 

Formule de Gauss 28 

Formule de Mellin 31 

Relation entre la fonction gamma et la série hypergéométrique 33 



IV. 



FONCTION DE BINET. 

Définition de la fonction de Binet 36 

Formules de Binet 38 

Formule de Gudermann 40 

Nouvelle expression du module de T(a + fii) 42 

Formule de Stirling 43 

Limites du nombre de Bernoulli B n 44 

Calcul d'une limite 45 



V. 



FONCTIONS O(x) ET W(x). 

Définition et propriétés des fonctions O(x) et W(x) 46 

/Il 1 \ 

Limite de l'expression log n — I - H + • • • H I pour n - oo 49 

\ JL JL \ L JL ~t~ fi L / 

Développement de la différence 0(x + a) - O(x) 49 

Détermination de O(x) pour une valeur rationnelle de x 53 

Étude de l'équation O(x) = 58 

Courbe figurative de la fonction gamma 60 



— 87 — 

VI. 

DÉVELOPPEMENTS EN SÉRIES ENTIÈRES. 



Développement de log T(l + x) en série entière 62 

1 

Développement de T(l + x) et de — en séries entières 65 

Développement de P(l + x) et de Q(l + x) en séries entières 67 

Développement de 0(1 + x) et de W(l + x) en séries entières 68 

VIL 



APPLICATIONS. 

Limite de produits infinis convergents dont le terme général est une fonction 

rationnelle de l'indice 70 

Sommation de séries dont le terme général est une fonction rationnelle de 

l'indice 73 

Résolution de certaines équations fonctionnelles 75 

Équation de Lindhagen 79 

Équations de Crelle 81 



88 



TABLE ALPHABETIQUE. 



Abel, 75 

Abel (Théorème d'), 56, 63 
Alexejewsky (Fonctions d'), 10 
Andoyer, 5, 21, 47, 73 
Appell, 4, 10, 16, 21, 70, 75 
Appell (Fonctions d'), 10, 70 
Appell (Problème d'), 21 

Barnes, 4, 9, 47 

Barnes (Fonction gamma double de), 9 

Bauer, 53 

Bellavitis, 65 

Berger, 7 

Bernoulli (Daniel), 1 

Bernoulli (Jean), 53 

Bernoulli (Limites des nombres de), 

44^5 
Bessel, v, 3, 4, 82 
Binet, 1-3, 28, 30, 37-39, 53 
Binet (Fonction de), 36-45 
Binet (Formules de), 38-40 
Binet (Identité de), 53 
Blaserna (Fonction de), 48 
Bourget, 22, 25, 38, 40, 59, 66-68 
Bourguet (Fonction de), 22 
Brunel, 4 

Cauchy, 3, 36, 45, 51 

Cayley, 44 

Cesàro, 7, 44 

Constante d'Euler, 8-9 

Constante d'Euler (Calcul de la), 8, 64 

Courbe figurative de la fonction 

gamma, 60-61 
Crelle, 4, 82 



Crelle (Équations de), 70, 81-84 

Développement de T(l + x) et de 
1 
— en séries entières, 10-11, 

r(i + x) 

65-67 
Développement de la différence 0(x + 

a) - O(x), 49-53 
Développement de log T(l +x) en série 

entière, 62-65 
Développement de P(l + x) et de Q(l + 

x) en séries entières, 21, 67-68 
Développement de 0(1 +x) et de ^(1 + 

x) en séries entières, 68 
Dirichlet, voir Lejeune-Dirichlet 

Enestrôm, 72 

Engelmann, 82 

Équation P(x) = (Étude de 1'), 23-25 

Équation O(x) = (Étude de 1'), 58-59 

Équations fonctionnelles (Résolution 

de certaines), 75-78 
Euler, v, 1, 2, 4, 26, 27, 30, 46, 63, 64 
Euler (Constante d'), 7-8, 64 
Euler (Constante d') , 49 
Eytelwein, 4 

Féaux, 39 

Facultés analytiques (Théorie des), 5, 

8, 82-84 
Fonction O(x) (Détermination pour 

une valeur rationnelle de x de la), 

53-58 
Fonction factorielle de Weierstrass, 4, 

8, 22, 32 



89 



90 



Fonction gamma (Courbe figurative 

de la, 60-61 
Fonction gamma (Décomposition en 

facteurs primaires de la), 9-10 
Fonction gamma (Définition de la), 6-8 
Fonction gamma (Développements en 

séries entières), 10-11, 62-65 
Fonction gamma (Historique de la), 1- 

5 
Fonction gamma double de Barnes, 9 
Fonction ô)(x), 36-45 
Fonction T(x), 22 
Fonctions inexplicables, 2 
Fonctions P(x) et Q(x), 5, 19-25, 67-68 
Fonctions O(x) et W(x), 46-61, 68, 73- 

75,78 
Fuss (P.-H.), 1 

Gasparis (A. de), 5, 72 

Gauss, v, vi, 2, 4, 7, 30, 34, 47, 51, 53, 

58, 64, 69 
Gauss (Formule de), 28-31, 45 
Gauss (Règle de convergence de), 73 
Genocchi, 40 
Gilbert, 40 
Glaisher, 8, 44 

Godefroy (Maurice), 11, 33, 67 
Goldbach, 1 
Gudermann, 4, 42 
Gudermann (Formule de), 40-42 

Hôlder, 47 

Hankel, v 

Heine (E.), 9 

Heine (Fonctions de), 9 

Hermite, 4, 5, 19, 21, 46, 47, 59, 73, 75 

Hoppe, 31 

Jeffery, 65 

Jensen, 4, 42, 53, 75, 79 

Kramp, 4, 82 
Kummer, 4 

Legendre, v, 3, 4, 7, 28, 31, 58, 63, 64 
Legendre (Formule de), 27-28, 30-31 



Lejeune-Dirichlet, 3 

Lerch (Mathias), 4, 13, 42, 53 

Lie (Sophus), 75 

Limite pour n - oo de logn — I — \- 
1 + ... + _L_\ 4 9 



x + 1 x+n—1 

Limite pour n - oo de 

1 F(x + n) M 

n x i-2---(n-l)' 
Lindhagen, 4, 22, 25, 67, 79 
Lindhagen (Équation de), 70, 79-81 
Liouville, v, vi, 3, 32 

Mansion, 44 

Mellin, v, 4, 5, 19, 21, 32, 34, 77 

Mellin (Fonction de), 21, 77 

Mellin (Formule de), 31-33 

Minimum de T(x), 59 

Mittag-Leffler (Théorème de), 5, 20, 21, 

46,73 
Module de T(x), 12-13, 42 

Newman, 9 
Nicolai, 64, 69 

Oettinger, 4, 82 
Ohm, 4, 82 

Picard (Emile), 4 

Pincherle (Fonction de), 10 

Plana, 37 

Prym, v, 4, 5, 14, 19 

Prym (Fonctions de), 5, 14, 17-25 

Pringsheim, v 

Produits infinis (Limite de), 70-73 

Résidus de T(x), 13-14 

Raabe, 3, 4 

Relation des compléments, 26-27, 33- 

35 
Relation entre la fonction gamma et la 

série hypergéométrique, 33-35 
Relation fonctionnelle de T(x), 11-12 
Rouché, 44 

Série hypergéométrique (Conver- 
gence de la), 16 



91 



Série hypergéométrique (Relation 
entre la fonction gamma et la), 
33-35 

Séries (Sommation de), 73-75 

Séries hypergéométriques d'ordre su- 
périeur, 77 

Scheefer, 20 

Scheibner, 70 

Schering, 3 

Shanks, 8 

Sonine, 30 

Stieltjes, 13, 63, 64 

Stirling, 1, 15, 16, 43, 72 

Stirling (Formule de), 42-45 

Stirling (Série de), 1 



Sylow, 75 

Tannery (Jules), 11 
Teixeira (Gomes), 44 
Thomae, 35 
Tilly (de), 40 

Vandermonde, 4, 82 

Wallis, 1 

Wallis (Formule de), 1, 28, 30, 41, 44 
Weierstrass, v, 4, 8-10, 14, 46, 82, 84 
Weierstrass (Fonction factorielle de), 4, 

8, 22, 32 
Weierstrass (Formule de), 8-10, 26 
Weyr (Eduard), 49, 75 



30254 



Paris. — Imprimerie GAUTHIER-VILLARS, quai des Grands-Augustins, 55. 



End of the Project Gutenberg EBook of La Fonction Gamma, by Maurice Godefroy 

*** END OF THIS PROJECT GUTENBERG EBOOK LA FONCTION GAMMA *** 

***** This file should be named 2980Q-pdf.pdf or 29800-pdf .zip ***** 
This and ail associated files of various formats will be found in: 
http : //www . gutenberg . org/2/9/8/0/29800/ 

Produced by Joshua Hutchinson, Andrew D. Hwang , and the 
Online Distributed Proofreading Team at http://www.pgdp.net 
(This file was produced from images from the Cornell 
University Library: Historical Mathematics Monographs 
collection.) 



Updated éditions will replace the previous one--the old éditions 
will be renamed. 

Creating the works from public domain print éditions means that no 
one owns a United States copyright in thèse works, so the Foundation 
(and you!) can copy and distribute it in the United States without 
permission and without paying copyright royalties. Spécial rules, 
set forth in the General Terms of Use part of this license, apply to 
copying and distributing Project Gutenberg-tm electronic works to 
protect the PROJECT GUTENBERG-tm concept and trademark. Project 
Gutenberg is a registered trademark, and may not be used if you 
charge for the eBooks, unless you receive spécifie permission. If you 
do not charge anything for copies of this eBook, complying with the 
rules is very easy. You may use this eBook for nearly any purpose 
such as création of derivative works, reports, performances and 
research. They may be modified and printed and given away--you may do 
practically ANYTHING with public domain eBooks. Redistribution is 
subject to the trademark license, especially commercial 
redistribution. 



*** START: FULL LICENSE *** 

THE FULL PROJECT GUTENBERG LICENSE 

PLEASE READ THIS BEFORE YOU DISTRIBUTE OR USE THIS WORK 

To protect the Project Gutenberg-tm mission of promoting the free 
distribution of electronic works, by using or distributing this work 
(or any other work associated in any way with the phrase "Project 
Gutenberg"), you agrée to comply with ail the terms of the Full Project 
Gutenberg-tm License (available with this file or online at 
http : //gutenberg . org/license) . 



Section 1. General Terms of Use and Redistributing Project Gutenberg-tm 
electronic works 

l.A. By reading or using any part of this Project Gutenberg-tm 
electronic work, you indicate that you hâve read, under stand, agrée to 
and accept ail the terms of this license and intellectual property 
(trademark/copyright) agreement. If you do not agrée to abide by ail 



the terms of this agreeraent, you must cease using and return or destroy 
ail copies of Project Gutenberg-tm electronic works in your possession. 
If you paid a fee for obtaining a copy of or access to a Project 
Gutenberg-tm electronic work and you do not agrée to be bound by the 
terras of this agreement, you may obtain a refund from the person or 
entity to whom you paid the fee as set forth in paragraph l.E.8. 

l.B. "Project Gutenberg" is a registered trademark. It may only be 
used on or associated in any way with an electronic work by people who 
agrée to be bound by the terms of this agreeraent. There are a few 
things that you can do with most Project Gutenberg-tm electronic works 
even without complying with the full terms of this agreement. See 
paragraph l.C below. There are a lot of things you can do with Project 
Gutenberg-tm electronic works if you follow the terms of this agreement 
and help préserve free future access to Project Gutenberg-tm electronic 
works. See paragraph l.E below. 

l.C. The Project Gutenberg Literary Archive Foundation ("the Foundation" 
or PGLAF) , owns a compilation copyright in the collection of Project 
Gutenberg-tm electronic works. Nearly ail the individual works in the 
collection are in the public domain in the United States. If an 
individual work is in the public domain in the United States and you are 
located in the United States, we do not claim a right to prevent you from 
copying, distributing , performing, displaying or creating derivative 
works based on the work as long as ail références to Project Gutenberg 
are removed. Of course, we hope that you will support the Project 
Gutenberg-tm mission of promoting free access to electronic works by 
freely sharing Project Gutenberg-tm works in compliance with the terms of 
this agreement for keeping the Project Gutenberg-tm name associated with 
the work. You can easily comply with the terras of this agreement by 
keeping this work in the same format with its attached full Project 
Gutenberg-tm License when you share it without charge with others. 

l.D. The copyright laws of the place where you are located also govern 
what you can do with this work. Copyright laws in most countries are in 
a constant state of change. If you are outside the United States, check 
the laws of your country in addition to the terms of this agreeraent 
before downloading, copying, displaying, performing, distributing or 
creating derivative works based on this work or any other Project 
Gutenberg-tm work. The Foundation makes no représentations concerning 
the copyright status of any work in any country outside the United 
States. 

l.E. Unless you hâve removed ail références to Project Gutenberg: 

l.E.l. The following sentence, with active links to, or other immédiate 
access to, the full Project Gutenberg-tm License must appear prominently 
whenever any copy of a Project Gutenberg-tm work (any work on which the 
phrase "Project Gutenberg" appears, or with which the phrase "Project 
Gutenberg" is associated) is accessed, displayed, performed, viewed, 
copied or distributed: 

This eBook is for the use of anyone anywhere at no cost and with 
alraost no restrictions whatsoever. You may copy it, give it away or 
re-use it under the terms of the Project Gutenberg License included 
with this eBook or online at www.gutenberg.org 



I.E.2. If an individual Project Gutenberg-tra electronic work is derived 
frora the public domain (does not contain a notice indicating that it is 
posted with permission of the copyright holder) , the work can be copied 
and distributed to anyone in the United States without paying any fées 
or charges. If you are redistributing or providing access to a work 
with the phrase "Project Gutenberg" associated with or appearing on the 
work, you must comply either with the requirements of paragraphs l.E.l 
through 1.E.7 or obtain permission for the use of the work and the 
Project Gutenberg-tm trademark as set forth in paragraphs 1.E.8 or 
l.E.9. 

I.E.3. If an individual Project Gutenberg-tm electronic work is posted 
with the permission of the copyright holder, your use and distribution 
must comply with both paragraphs l.E.l through 1.E.7 and any additional 
terms imposed by the copyright holder. Additional terms will be linked 
to the Project Gutenberg-tm License for ail works posted with the 
permission of the copyright holder found at the beginning of this work. 

I.E.4. Do not unlink or detach or remove the full Project Gutenberg-tm 
License terms from this work, or any files containing a part of this 
work or any other work associated with Project Gutenberg-tm. 

I.E.5. Do not copy, display, perform, distribute or redistribute this 
electronic work, or any part of this electronic work, without 
prominently displaying the sentence set forth in paragraph l.E.l with 
active links or immédiate access to the full terms of the Project 
Gutenberg-tm License. 

I.E.6. You may convert to and distribute this work in any binary, 
compressed, marked up, nonproprietary or proprietary form, including any 
word processing or hypertext form. However, if you provide access to or 
distribute copies of a Project Gutenberg-tm work in a format other than 
"Plain Vanilla ASCII" or other format used in the officiai version 
posted on the officiai Project Gutenberg-tm web site (www.gutenberg.org), 
you must, at no additional cost, fee or expense to the user, provide a 
copy, a means of exporting a copy, or a means of obtaining a copy upon 
request, of the work in its original "Plain Vanilla ASCII" or other 
form. Any alternate format must include the full Project Gutenberg-tm 
License as specified in paragraph l.E.l. 

I.E.7. Do not charge a fee for access to, viewing, displaying, 
performing, copying or distributing any Project Gutenberg-tm works 
unless you comply with paragraph 1.E.8 or l.E.9. 

I.E.8. You may charge a reasonable fee for copies of or providing 
access to or distributing Project Gutenberg-tm electronic works provided 
that 

- You pay a royalty fee of 28% of the gross profits you dérive from 

the use of Project Gutenberg-tm works calculated using the method 
you already use to calculate your applicable taxes. The fee is 
owed to the owner of the Project Gutenberg-tm trademark, but he 
has agreed to donate royalties under this paragraph to the 
Project Gutenberg Literary Archive Foundation. Royalty payments 
must be paid within 6Q days following each date on which you 
prépare (or are legally required to prépare) your periodic tax 
returns . Royalty payments should be clearly marked as such and 



sent to the Project Gutenberg Literary Archive Foundation at the 
address specified in Section 4, "Information about donations to 
the Project Gutenberg Literary Archive Foundation." 

- You provide a full refund of any money paid by a user who notifies 

you in writing (or by e-mail) within 30 days of receipt that s/he 
does not agrée to the terms of the full Project Gutenberg-tm 
License . You must require such a user to return or 
destroy ail copies of the works possessed in a physical médium 
and discontinue ail use of and ail access to other copies of 
Project Gutenberg-tm works. 

- You provide, in accordance with paragraph 1.F.3, a full refund of any 

money paid for a work or a replacement copy, if a defect in the 
electronic work is discovered and reported to you within 90 days 
of receipt of the work. 

- You comply with ail other terms of this agreement for free 

distribution of Project Gutenberg-tm works. 

I.E.9. If you wish to charge a fee or distribute a Project Gutenberg-tm 
electronic work or group of works on différent terms than are set 
forth in this agreement, you must obtain permission in writing from 
both the Project Gutenberg Literary Archive Foundation and Michael 
Hart, the owner of the Project Gutenberg-tm trademark. Contact the 
Foundation as set forth in Section 3 below. 

l.F. 

l.F.l. Project Gutenberg volunteers and employées expend considérable 
effort to identify, do copyright research on, transcribe and proofread 
public domain works in creating the Project Gutenberg-tm 
collection. Despite thèse efforts, Project Gutenberg-tm electronic 
works, and the médium on which they may be stored, may contain 
"Defects," such as, but not limited to, incomplète, inaccurate or 
corrupt data, transcription errors, a copyright or other intellectual 
property infringement , a defective or damaged disk or other médium, a 
computer virus , or computer codes that damage or cannot be read by 
your equipment . 

l.F. 2. LIMITED WARRANTY, DISCLAIMER OF DAMAGES - Except for the "Right 
of Replacement or Refund" described in paragraph l.F. 3, the Project 
Gutenberg Literary Archive Foundation, the owner of the Project 
Gutenberg-tm trademark, and any other party distributing a Project 
Gutenberg-tm electronic work under this agreement, disclaim ail 
liability to you for damages, costs and expenses, including légal 
fées. YOU AGREE THAT YOU HAVE NO REMEDIES FOR NEGLIGENCE, STRICT 
LIABILITY, BREACH OF WARRANTY OR BREACH OF CONTRACT EXCEPT THOSE 
PROVIDED IN PARAGRAPH F3. YOU AGREE THAT THE FOUNDATION, THE 
TRADEMARK OWNER, AND ANY DISTRIBUTOR UNDER THIS AGREEMENT WILL NOT BE 
LIABLE TO YOU FOR ACTUAL, DIRECT, INDIRECT, CONSEQUENTIAL , PUNITIVE OR 
INCIDENTAL DAMAGES EVEN IF YOU GIVE NOTICE OF THE POSSIBILITY OF SUCH 
DAMAGE . 

l.F. 3. LIMITED RIGHT OF REPLACEMENT OR REFUND - If you discover a 
defect in this electronic work within 90 days of receiving it, you can 
receive a refund of the money (if any) you paid for it by sending a 



written explanation to the person you received the work from. If you 
received the work on a physical médium, you must return the médium with 
your written explanation. The person or entity that provided you with 
the defective work may elect to provide a replacement copy in lieu of a 
refund. If you received the work electronically , the person or entity 
providing it to you may choose to give you a second opportunity to 
receive the work electronically in lieu of a refund. If the second copy 
is also defective, you may demand a refund in writing without further 
opportunities to fix the problem. 

I.F.4. Except for the limited right of replacement or refund set forth 
in paragraph 1.F.3, this work is provided to you 'AS-IS' WITH NO OTHER 
WARRANTIES OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO 
WARRANTIES OF MERCHANTIBILITY OR FITNESS FOR ANY PURPOSE. 

I.F.5. Some states do not allow disclaimers of certain implied 
warranties or the exclusion or limitation of certain types of damages. 
If any disclaimer or limitation set forth in this agreement violâtes the 
law of the state applicable to this agreement, the agreement shall be 
interpreted to make the maximum disclaimer or limitation permitted by 
the applicable state law. The invalidity or unenforceability of any 
provision of this agreement shall not void the remaining provisions. 

I.F.6. INDEMNITY - You agrée to indemnify and hold the Foundation, the 
trademark owner, any agent or employée of the Foundation, anyone 
providing copies of Project Gutenberg-tm electronic works in accordance 
with this agreement, and any volunteers associated with the production, 
promotion and distribution of Project Gutenberg-tm electronic works, 
harmless from ail liability, costs and expenses, including légal fées, 
that arise directly or indirectly from any of the following which you do 
or cause to occur: (a) distribution of this or any Project Gutenberg-tm 
work, (b) altération, modification, or additions or deletions to any 
Project Gutenberg-tm work, and (c) any Defect you cause. 



Section 2. Information about the Mission of Project Gutenberg-tm 

Project Gutenberg-tm is synonymous with the free distribution of 
electronic works in formats readable by the widest variety of computers 
including obsolète, old, middle-aged and new computers. It exists 
because of the efforts of hundreds of volunteers and donations from 
people in ail walks of life. 

Volunteers and financial support to provide volunteers with the 
assistance they need, are critical to reaching Project Gutenberg-tm' s 
goals and ensuring that the Project Gutenberg-tm collection will 
remain freely available for générations to corne. In 2001, the Project 
Gutenberg Literary Archive Foundation was created to provide a secure 
and permanent future for Project Gutenberg-tm and future générations. 
To learn more about the Project Gutenberg Literary Archive Foundation 
and how your efforts and donations can help, see Sections 3 and 4 
and the Foundation web page at http://www.pglaf.org. 



Section 3. Information about the Project Gutenberg Literary Archive 
Foundation 



The Project Gutenberg Literary Archive Foundation is a non profit 
5Ql(c)(3) educational corporation organized under the laws of the 
state of Mississippi and granted tax exempt status by the Internai 
Revenue Service. The Foundation' s EIN or fédéral tax identification 
nuraber is 64-6221541. Its 501(c)(3) letter is posted at 
http://pglaf.org/fundraising. Contributions to the Project Gutenberg 
Literary Archive Foundation are tax déductible to the full extent 
perraitted by U.S. fédéral laws and your state 's laws. 

The Foundation' s principal office is located at 4557 Melan Dr. S. 
Fairbanks, AK, 99712., but its volunteers and employées are scattered 
throughout numerous locations. Its business office is located at 
809 North 1500 West, Sait Lake City, UT 84116, (801) 596-1887, email 
business@pglaf.org. Email contact links and up to date contact 
information can be found at the Foundation' s web site and officiai 
page at http://pglaf.org 

For additional contact information: 
Dr. Gregory B. Newby 
Chief Executive and Director 
gbnewby@pglaf.org 



Section 4. Information about Donations to the Project Gutenberg 
Literary Archive Foundation 

Project Gutenberg-tm dépends upon and cannot survive without wide 
spread public support and donations to carry out its mission of 
increasing the number of public domain and licensed works that can be 
freely distributed in machine readable form accessible by the widest 
array of equipment including outdated equipment . Many small donations 
($1 to $5,000) are particularly important to maintaining tax exempt 
status with the 1RS. 

The Foundation is committed to complying with the laws regulating 
charities and charitable donations in ail 58 states of the United 
States. Compliance requirements are not uniform and it takes a 
considérable effort, much paperwork and many fées to meet and keep up 
with thèse requirements. We do not solicit donations in locations 
where we hâve not received written confirmation of compliance. To 
SEND DONATIONS or détermine the status of compliance for any 
particular state visit http://pglaf.org 

While we cannot and do not solicit contributions from states where we 
hâve not met the solicitation requirements, we know of no prohibition 
against accepting unsolicited donations from donors in such states who 
approach us with offers to donate. 

International donations are gratefully accepted, but we cannot make 
any statements concerning tax treatment of donations received from 
outside the United States. U.S. laws alone swamp our small staff. 

Please check the Project Gutenberg Web pages for current donation 
methods and addresses. Donations are accepted in a number of other 
ways including checks, online payments and crédit card donations. 
To donate, please visit: http://pglaf.org/donate 



Section 5. General Information About Project Gutenberg-tm electronic 
works . 

Professor Michael S. Hart is the originator of the Project Gutenberg-tm 
concept of a library of electronic works that could be freely shared 
with anyone. For thirty years, he produced and distributed Project 
Gutenberg-tm eBooks with only a loose network of volunteer support. 



Project Gutenberg-tm eBooks are often created from several printed 
éditions, ail of which are confirmed as Public Domain in the U.S. 
unless a copyright notice is included. Thus, we do not necessarily 
keep eBooks in compliance with any particular paper édition. 



Most people start at our Web site which has the main PG search facility: 

http : //www . gutenberg . org 

This Web site includes information about Project Gutenberg-tm, 
including how to make donations to the Project Gutenberg Literary 
Archive Foundation, how to help produce our new eBooks, and how to 
subscribe to our email newsletter to hear about new eBooks. 



